Electron transport in disordered Luttinger liquid 
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We study the transport properties of interacting electrons in a disordered quantum wire within 
the framework of the Luttinger liquid model. We demonstrate that the notion of weak locaUzation is 
applicable to the strongly correlated one-dimensional electron system. Two alternative approaches 
to the problem are developed, both combining fermionic and bosonic treatment of the underlying 
physics. We calculate the relevant dephasing rate, which for spinless electrons is governed by 
the interplay of electron-electron interaction and disorder, thus vanishing in the clean limit. Our 
approach provides a framework for a systematic study of mesoscopic effects in strongly correlated 
electron systems. 

PACS numbers: 71.10.Pm, 73.21.-b, 73.63.-b, 73.20.Jc 



I. INTRODUCTION 

Mesoscopics of strongly correlated electron systems 
has emerged as an area of great interest to both ex- 
perimental and theoretical communities working in the 
field of nanoscale physics. Recently, progress in man- 
ufacturing of nanodevices has paved the way for sys- 
tematic transport measurements on narrow quantum 
wires with a few or single conducting channels. Most 
prominent examples of these are: single-wall carbon 
nanotubeSfi. cleaved-edge^, TZ-groovCf^ and crystallized 
in a matrisi semiconductor quantum wires, quantum 
Hall edges running in opposite directions and intercon- 
nected by means of tunneling, polymer nanofiberSji 
and metallic nanowires.^^ 

Much experimental attention has been focused on the 
electrica l prop erties of single- wall nanotubes (RefsH^III!, 
[illilliil^IiBli^Iiililafllklillllll^ and references 
therein). In particular, evidence has emerged pointing to- 
wards the existence of Luttinger liquid in metallic single- 
wall carbon nanotubes, as expected for strongly inter- 
acting electrons in one dimension (Id). The Luttinger 
liquid behavior was observed via the power-law tempera- 
ture and bias- voltage dependence of the current through 
tunneling contacts attached to the nanotubes. In the 
past few years, technological advances have made pos- 
sible the fabrication of contacts between nanotubes and 
metallic leads with a very low contact resistance, which 
allows one to observe mesoscopic effects J^i^SiSii^ Fur- 
ther, techniques to grow very long (of the mm scale) 
nanotubes have been developed and the corresponding 
transport measurements performedi^ 

On the theoretical side, the challenge is to expand 
the ideas that have been developed for mesoscopic dis- 
ordered systems on one side and for strongly correlated 
clean systems on the other. The situation is particu- 
larly interesting in Id, where both disorder and electron- 
electron (e-e) interaction, even if they are weak, modify 
dramatically the large-scale, low-energy physics of the 
problem. In a clean wire, e-e interaction leads to the 
formation of a Luttinger-liquid ground state. The prop- 



erties of the Luttinger liquid without impurities and in 
the presence of a single barrier are known in great de- 
tail, see Refs. for review. Much less 
is known about the Luttinger liquid in the presence of 
many impurities (disordered Luttinger liquids). The in- 
fluential works by Apel and Rice^^ and by Giamarchi and 
SchulaSi have defined the state of the art in this area for 
two decades. Recent years have seen a revival of inter- 
est in disordered Luttinger liquids ^i^i^ii^i^adaiild^d^ 
largely motivated by the technological advances men- 
tioned above. However, the very applicability of such 
key notions of the mesoscopic physics as weak localization 
(WL) and dephasing to a disordered Id system, not to 
mention a detailed analysis of these effects, remained an 
open problem. This problem is a subject of the present 
paper. Throughout the paper we study spinless (spin- 
polarized) electrons. Spin- related effects will be consid- 
ered in Ref. [ij . 

The structure of the paper is as follows. In Sec. |TT] 
we summarize the earlier achievements in the mesoscop- 
ics of higher-dimensional systems on one side and in the 
theory of strongly correlated Id systems on the other. 
Section IIIII contains an exact formulation of the prob- 
lem, as well as a discussion of our strategy in solving 
it. In Sec. IIVI we present the bosonization method and 
the renormalization-group technique that we use to ac- 
count for the Luttinger-liquid renormalizations. Having 
analyzed in Sec. |V] the e-c inelastic scattering rate in Id, 
we turn to the calculation of a weak-localization correc- 
tion to the conductivity and the associated dephasing 
rate in Sec. |Vll We also compare this dephasing rate 
with the one that governs the damping of Aharonov- 
Bohm oscillations. Section IVIII is devoted to the anal- 
ysis of the transport properties of a disorder Luttinger 
liquid by using an alternative approach — the "functional 
bosonization" — which allows one to treat the renormal- 
ization and the inelastic scattering on an equal footing. 
Our results are summarized in Sec. IVIIII Technical de- 
tails of the calculations are presented in several Appen- 
dices. Some of the results of this paper were published 
in a brief form in Letter Ref. 
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II. BACKGROUND 

In this introductory section, we briefly summarize the 
known results for higher-dimensional mesoscopic sys- 
tems (Sec. Ill A|) and for strongly correlated Id systems 
(Sec. IIIBI) . These will serve as a starting point for our 
work. 



A. Mesoscopics of higher-dimensional systems 

During the last three decades, mesoscopic physics of 
low-dimensional diffusive systems has been attracting a 
great deal of attention . It has been recog- 
nized that the effects of the quantum interference and e-e 
interaction become strongly enhanced in disordered sys- 
tems of reduced dimensionality d, affecting in an essen- 
tial way the low-temperature transport properties of the 
systems. These effects lead to a set of remarkable phe- 
nomena, including the WL, mesoscopic conductance fluc- 
tuations, and the Altshuler-Aronov interaction-induced 
quantum correction to the conductivity. At still lower 
temperatures, these effects may drive the system towards 
strong localization. We now briefly remind the reader the 
basic facts in this field. 



1. Moderately low temperatures: Weak-localization regime 

The quantum localization is the most prominent man- 
ifestation of the quantum interference in disordered sys- 
tems. At not too low temperatures, the localization ef- 
fects are cut off by inelastic processes, resulting in a WL 
correction^ Aa^ji to the Drude conductivity ctd, 
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where D is the diffusion constant, the density of states, 
I the mean free path, r the transport mean free time, 
T0 and = [Dt^)^/'^ are the interaction-induced de- 
phasing time and length, respectively. We use notation 
(dq) = dfiq/{2nY. Throughout the paper h—1. Another 
manifestation of the quantum coherence are mesoscopic 
fluctuations of the conductance and of other observables. 

As far as e-e interaction in the WL regime is concerned, 
it is responsible for two distinctly different main effects. 
First, it renormalizes elastic scattering on impurities 
through the creation of virtual electron-hole excitations 
(screening) with a high energy transfer lo larger than the 
temperature T, T < \uj\ < ep. Here ep is the Fermi 
energy. At sufficiently high T (in the ballistic regime 
Tt ^ 1), i.e., when the relevant spatial scales are smaller 
than Z, the screening can be described in terms of a dress- 
ing of impurities by Friedel oscillations^i^ and yields a 



T-dependent Drude conductivity.^^ At lower T (in the 
diffusive regime Tt ^ 1), effects of this type generate 
the Altshuler-Aronov correction to the conductivity,— 
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where Lt (D/TY^'^ is the thermal length. This cor- 
rection is perturbative both in the effective interaction 
constant a and in the disorder strength. The Coulomb 
interaction corresponds to a ~ 1. The T dependence 
of AdAA is determined by the thermal smearing of the 
distribution function, which introduces an infrared cut- 
off on spatial scales of the order of Lt- Using the 
renormalization-group (RG) approach within the repli- 
cated tJ-model, Finkel'stein extended this theory to in- 
clude the interaction nonperturbatively;^ 

Effects related to the interaction-induced renormaliza- 
tions and screening are also responsible for a zero-bias 
anomaly^i^ in the tunneling density of states as a func- 
tion of electron energy, i'{e). It is worth noting that, 
in the case of Coulomb interaction, ^{e) may become 
strongly suppressed , ^^i^^i^^ 
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where k^^ is the screening radius, while the conductivity 
can still remain close to its Drude value (Actaa/o'd <C 1). 

The second effect of interaction is inelastic e-e scatter- 
ing which breaks the phase coherence.— A key concept 
in the localization theory of a disordered Fermi liquid 
is that of the dephasing rate r^^. The phase-breaking 
processes are of crucial importance in the problem of 
transport since without dephasing the conductivity of 
low-dimensional systems would be zero at any T, see 
Sec. Ill A^ Two qualitatively different sources of dephas- 
ing are possible: (i) scattering of electrons by external 
excitations (in practice, phonons) and (ii) e-e scattering 
(at low T the dephasing is mostly due to e-e interactions). 
In either case, at sufficiently high T, the dephasing rate 
is high, so that the localization effects are reduced 
to the WL correction to the Drude conductivity. 

A characteristic energy transfer in the real inelas- 
tic processes that determine the phase relaxation is re- 
stricted by temperature: \lu\ < T, so that the dephasing 
rate vanishes at zero temperature. The infrared diver- 
gency characteristic of a diffusive system makes it neces- 
sary to introduce a self-consistent cutoff in the equation 
for the dephasing rate^— 
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A rigorous calculation necessitates employing the path- 
integral techniquei^ 

The different energy scales relevant to the disorder 
renormalization due to virtual processes with T < |w| < 
ep and to the dephasing due to real inelastic scattering 
with ^ T allow for a straightforward separation of 
these two effects of e-e interaction. At sufficiently high 
T, the conductivity 



ct(T) ~ (7d + AcTwl + AfJAA 



(2.5) 



is close to the Drude value ctd- It is worth mentioning 
that for weak interaction (a ^ 1) the WL correction 
[Eq. (|2.ip with given by Eq. (|2.4p ] is much stronger 
than the Altshuler-Aronov correction [Eq. (|2.2p ]. The 
WL correction grows with lowering T and eventually be- 
comes strong (|A(Twi|/o'D ~ 1) when reaches ^ or, 
equivalently, when becomes of order A^. Here ^ is 
the localization length and A^ is the corresponding en- 
ergy level spacing. At this temperature, IActaaI/cd ^ 1 
for a <C 1. This means that in a weakly interacting sys- 
tem the strong localization occurs due to the growth of 
the interference-induced WL correction. In the next sub- 
section, we briefly overview the strong localization regime 
which emerges at low temperatures, when t7^ < A^. 



2. Low temperatures: Fate of dephasing in the strongly 
localized regime 

In a pathbreaking paper— Anderson demonstrated 
that a quantum particle may become localized by a ran- 
dom potential. In particular, in noninteracting systems 
of one- or two-dimensional geometry even weak disorder 
localizes all electronic states,— thus leading to the ex- 
actly zero conductivity, a(T) = 0, whatever temperature 
T. A nonzero cr(T) in such systems may only occur due 
to inelastic scattering processes leading to the dephasing 
of electrons. 

As discussed in Sec. Ill A 1[ the WL correction p.ip 
to the Drude conductivity of a diffusive system behaves 



(2-d)/2 
^0 



(In T(j, for d ^ 2) and diverges with lowering 
temperature for d < 2. This singular behavior of the 
perturbative corrections is a manifestation of strong An- 
derson localization. It prompts a question as to how the 
disordered interacting low-dimensional systems conduct 
at low temperatures, when the localization effects are no 
longer weak. 

It is commonly believed that transport in the local- 
ized regime is of the hopping nature. Mott^ has de- 
veloped the picture of variable-range-hopping (VRH) as 
the transport mechanism in the insulating state. The 



VRH mechanism of Ref. lei" is well justified in the pres- 
ence of delocalized phonons which provide inelastic scat- 
tering between single-particle states with different ener- 
gies (phonon-assisted hopping). The finite-temperature 
conductivity is then proportional to the strength of 
the electron-phonon coupling and governed by Mott's 
VRH:6i ^(y) ^ exp[-(To/T)^] with ^ = l/{d+l) de- 
pending on the spatial dimensionality d. In the presence 
of a long-range Coulomb interaction, the Coulomb gap 
in the tunneling density of states modifies^ the VRH 
exponent, /i = i. 

But what if the coupling to phonons is negligibly weak? 
It was proposed in Ref. l63| that the e-e interaction by it- 
self is suflacient to induce VRH at low T. Although the 
idea^^ of a phononless VRH was widely exploited in the 
community, it is, however, in conflict with the argument 
by Fleishman and Anderson^ that elementary hops in 
the low-T limit are forbidden for a short-range interac- 
tion (including the case of 1/r Coulomb interaction in 
d < 3), because energy conservation cannot be respected 
when an electron attempts a real transition by exciting an 
electron-hole pair. Since there is a mismatch in the posi- 
tions of the energy levels of individual single-particle lo- 
calized states, inelastic scattering is necessary to broaden 
these levels to allow for hopping between them. It was 
argued in Ref. '6^ that at low T a short-range e-e interac- 
tion alone is not sufficient to broaden the localized states 
and to induce a nonzero conductivity. The situation is 
particularly interesting in Id and 2d, where no mobil- 
ity edge exists, activation to which otherwise might give 
a{T) ^0. 

Very recently, the problem was reconsidered in 
Rcfs. 65,66. It was found there that the Anderson lo- 
calization in Fock space crucially affects the conductivity 
at low T. As a result, the phononless VRH should not ex- 
ist; instead, the system undergoes a transition at certain 
temperature Tc below which the conductivity vanishes, 
in agreement with Ref. [6^ 

The mechanism of low-T transport in low-dimensional 
disordered systems is thus completely determined by in- 
elastic processes which govern the lifetime of localized 
states. This lifetime may be viewed as a low-temperature 
reincarnation of the dephasing time that controls the cut- 
off of quantum-interference effects at higher T. There- 
fore, an analysis of the WL effects turns out to be impor- 
tant also for understanding of the establishment of the 
strong localization with lowering T. 



B. Strongly correlated Id electron systems 

This paper is concerned with transport in a single- 
channel quantum wire, where both disorder and inter- 
action more strongly modify the transport properties as 
compared to higher dimensions. In the absence of disor- 
der, the e-e correlations drive the system into the non- 
Fermi liquid state known as Luttinger liquid. Yet another 
peculiarity of the single-channel Id system is that the bal- 
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listic motion on short scales crosses over in the absence 
of interaction directly to the localization regime, with no 
diffusive dynamics on intermediate scales. We are now 
going to expand on these known features of Id electronic 
systems. 

1. Noninteracting case: Anderson Localization 

We first briefly recall the peculiarities of Id noninter- 
acting disordered systems. In the absence of interaction, 
transport in a disordered single-channel Id gas has been 
studied in great detail,S>S'^'2ai2i>2^'23^'2^ It is known 
that arbitrarily weak disorder leads to the localization of 
all electronic states in this system, with the localization 
length ^ being the mean free path. The dc conductivity 
is then zero, cr(T) = for any T, since the temperature 
merely determines the distribution function over the lo- 
calized states. 

The characteristic spatial and temporal scales of the 
problem are probed in the absence of interactions by the 
ac response at external frequency Q. The appropriate 
technique for the calculation of the response, which ex- 
ploits the separation of fast and slow degrees of freedom 
in the case of weak disorder, epT 3> 1, was put forward 
in a pioneering paper by Berezinskii;®^ later, alternative 
approaches were developed by Abrikosov and Ryzhkin^^ 
and by Berezinskii and Gor'kov. ''■^ In the limit of low fre- 
quency, Qt <^ 1, the leading contribution to the ac con- 
ductivity is given by the Berezinskii-Mott formulai ^^i^° 

Rea{fl)/aB ^ H^rfin^i^T) . (2.6) 

In the opposite (ballistic) limit of high frequency, fir ^ 
1, the conductivity a{Q) is approximately equal to the 
Drude formula, auiS^) = o'd/(1 — ifir), with a quantum 
interference (WL) correction arising at fourth order in 
the expansion in (17r)~^ <C 1. At fir ~ 1, the ballis- 
tic regime crosses over directly to the strongly localized 
regime, Eq. (|2.6p . so that the diffusive regime is absent in 
the single-channel quantum wire without e-e interaction. 

2. Clean quantum wires with interaction: Luttinger liquid 

The interest in the role of e-e interaction in quantum 
wires is largely inspired by the fact that in Id the inter- 
action breaks down the Landau's Fermi liquid approach. 
Owing to the particular geometry of the Fermi surface, 
systems of dimensionality one are unique in that the 
e-e correlations change the noninteracting picture com- 
pletely and lead to the formation of a strongly correlated 
state. A remarkable example of the correlated Id elec- 
tron phase is the Luttinger-liquid model (for a review 
see, e.g., Refs. [?7ll28ll2l30ll3l|ml). The Luttinger-liquid 
correlations show up in a power-law singularity of the 
tunneling density of states at low energy e counted from 
the Fermi surface, 

Ke)~^o(kl/A)'' . (2.7) 



Here Uq is the density of states in the absence of interac- 
tion, A is the ultraviolet cutoff (determined by the band 
width, the Fermi energy, or the range of interaction). 
The exponent /3 is determined by the interaction strength 
and depends on the system geometry (tunneling into the 
bulk of a wire differs from tunneling into the end of a 
wire) . The suppression of the tunneling density of states 
in Luttinger liquid is similar to the interaction-induced 
zero-bias anomaly in higher dimensions, Eq. (|2.3p . 

According to a common wisdom, electronic excitations 
do not exist in Luttinger liquid and a proper language 
to describe the strongly correlated electrons is the lan- 
guage of bosonic excitations. This notion is based on the 
fact that the bosonic density-fields in a clean Luttinger 
liquid do not decay in the linearized-spectrum approxi- 
mation, when a finite curvature of the electron spectrum 
is neglected. 



3. Renormalization of an impurity in Luttinger liquid 

Interactions between oppositely moving electrons gen- 
erate charge- and spin-density wave correlations that 
lead to striking transport properties of Luttinger liq- 
uid in the presence of impurities. The density- 
wave correlations yield a strong renormalization of an 
impurity; ''^I'^'^'i'^'^i'^^'^^i^"'^ if the bare strength of the impu- 
rity is small, the renormalized backscattering amplitude 
U{2kF, e) seen by an electron with energy e scales as 

U{2kF,e)^iA/\e\r , (2.8) 

where a is determined by the interaction strength. As a 
result, a single impurity in the infinite Luttinger liquid 
with a repulsive interaction effectively decouples the wire 
into two independent pieces, thus leading to a vanishing 
linear conductance at zero temperature. This conclusion 
remains valid no matter how weak the interaction is. 



4- Many impurities: Disordered Luttinger liquid 

The conductivity of a disordered Luttinger liquid was 
discussed for the first time more than 30 years ago in 
Refs. [tgIIttIItsI in terms of the renormalization of individ- 
ual impurities (see also Refs. Is^lssl for early attempts to 
construct the phase diagram of a disordered Luttinger 
liquid). In the beginning of the 80's, Apel and Rice^ 
combined the scaling arguments developed for noninter- 
acting systems^!^!^ with the results of the theory of a 
clean Luttinger liquid. Later, the problem of transport 
in a disordered Luttinger liquid was studied by perturba- 
tive in disorder methods based on bosonic field theories 
in Refs. I34|j83 l86l. A small correction to the quantized 
conductance of a short quantum wire due to a Gaussian 
disorder was derived in Ref. (for a multichannel wire 
see Ref. [st!) . A detailed exposition of these issues can be 
found in the recent lecture course Ref. [H. 



5 



The RG equations for a disordered Luttinger liquid 
were derived by Giamarchi and Schulz."^^ Their approach 
bears a similarity with the Finkel'stein RG^ developed 
earlier for higher-dimensional diffusive systems (both RG 
schemes treat disorder perturbatively and e-e interaction 
exactly). The flow equations of Ref. 113 describe, in ad- 
dition to the renornialization of impurities by interaction 
fSec. IIIlTSl) . also a renornialization of the interaction by 
disorder. A sufficiently strong attractive interaction was 
shown to yield a metal-insulator transitional at zero T 
with changing interaction strength. A repulsive interac- 
tion was predicted to enhance the disorder-induced lo- 
calization. In the latter case, the RG equations show 
that at a certain temperature (such that Lt ^ the 
renormalized disorder becomes strong. This means that 
the disordered interacting Id system enters a "strong- 
coupling regime" , where the Drude-like approach breaks 
down. At higher temperatures, these equations describe 
the temperature dependence of an effective scattering 
time renormalized by interaction (which can be viewed 
as screening of impurities by Friedel oscillations), thus 
yielding a temperature-dependent Drude conductivity, 

aD(T)oc(T/A)'". (2.9) 

Recently, a crossover between Eqs. (j2.9|) and (|2.2p with 
increasing number of channels was studied in Ref. 43. 

It is important to emphasize that the RG approach 
formalizes the renormalization of the Drude conductivity 
due to the T-dependent screening but does not capture 
the essential physics of the Anderson localization. In 
particular, the RG equations of Ref. 33 miss the inter- 
ference effects (coherent scattering on several impurities) 
that lead to the WL correction to the conductivity. As 
we demonstrate below, it is the WL correction that gov- 
erns the temperature dependence of the conductivity at 
moderately high temperatures for the case of weak inter- 
action and eventually drives the quantum wire into the 
strong-localization regime with decreasing temperature. 
Another related deficiency of the RG scheme of Ref. 's^ 
is the treatment of the scalings with length and u/T as 
interchangeable. While this approach is justified for the 
"elastic renormalization" of disorder (i.e., the screening 
effects), it does not properly account for the dephasing 
(since the dephasing length is in general different from 
Lt, see Sec. |IIIBlJand Sec UVCj) . 



5. Dephasing and localization in disordered Luttinger liquid 

As already mentioned in Sec. lIIB"4l a key ingredient of 
transport theory as regards the WL and the interaction- 
induced dephasing in a strongly correlated Id system is 
missing. Some papers (e.g., Refs. Is^lsih have suggested 
that the dephasing length that controls localization ef- 
fects in a disordered Luttinger liquid is Lt — u/T ^ where 
u is the plasmon velocity. According to this approach, the 
interference effects get strong with lowering T at Lt ~ 



An alternative approaches, is predicated on the assump- 
tion that the dephasing rate is determined by the single- 
particle properties of a clean Luttinger liquid. On top 
of that, one might think that since in the case of lin- 
ear dispersion the clean interacting electron system can 
be equivalently represented in terms of noninteracting 
bosons, the interaction should not induce any dephasing 
at all. The conductivity would then be exactly zero at 
any T. As we demonstrate below, none of the approaches 
captures the essential physics of dephasing in the conduc- 
tivity of a disordered Id system. 

As for the strongly localized regime, recently Ref. 113 
used bosonization to study the problem and concluded 
that transport is of variable-range-hopping character. 
Subsequently, Ref. ^ argued that this result requires a 
nonzero coupling to a thermal bath, in agreement with 
Ref. m, see Sec. HTXl 



C. Summary of known results: Id vs 2d 

We summarize the basic facts discussed in Sees. Ill Al 
and III Bl in Table |T| for simplicity, we compare here 2d 
and single-channel Id systems. It turns out that the dif- 
ference between the disordered Luttinger liquid and the 
higher-dimensional (quasi-ld and 2d) systems is not so 
strong as one might think judging by the proverbial non- 
Fermi-liquid character of the Luttinger liquid. 



III. DISORDERED LUTTINGER LIQUID 
A. Formulation of the problem 

Let us now introduce the model we are going to study. 
We consider a single-channel infinite wire at sufficiently 
low temperatures, much lower than the bandwidth and 
the Fermi energy. The effective low-energy theory can 
then be described by the Luttinger model )^''i^^'^^i'^°'^^i'^^ 
with a linear dispersion of the electronic spectrum Ck = 
kvp, where vp is the Fermi velocity. We use the notion of 
right- and left- movers corresponding to two Fermi points 
k ^ -iikp. The fermionic operators V'/i are then labeled 
hy fjL — ± for the right-/left-movers. 

We consider a finite-range (screened) pairwise e-e inter- 
action potential V{x—x'). We need only its Fourier trans- 
forms at g = and q = ±2kp: Vf and Vb, respectively, 
which give the amplitudes of forward and backward scat- 
tering. In terms of the g-ologyj^^i^SiiH the interaction in 
a homogeneous system is characterized by three coupling 
constants gi — Vb (e-e backscattering) , 172 = Vf (forward 
scattering between right- and left-movers), and 174 = Vf 
(forward scattering of right-/left-movers on right-/left- 
movers). In several places below, it will be convenient to 
treat g2 and 54 as independent parameters. 

Next, we introduce weak {epTo 3> 1) white-noise dis- 
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TABLE I: The 2d- Id vocabulary. 


2d 


Id 


disorder, no interaction: 


disorder, no interaction: 


cr(J ) - U 


a(i ) - U 


interaction, common tool: 


interaction, common tool: 


lermionic diagrams 


bosonization 


field theory: a-model 


field theory: 


Q-matrices 


bosonic density fields 


zero-bias anomaly in 


Luttinger ("non-Fermi") liquid, 


tunnelmg Dob, bq. (|2.3p 


power- law Dob, Eq. (|2.7|l 


T-dependent screening 


renormalization of impurity 


(Friedel oscillations) 


(Friedel oscillations), Eq. (12. 8p 


Altshuler-Aronov correction 


renormalization of disorder 


to the conductivity, Eq. (|2.2|l 


— * au(T), Eq. (fO)) 


Finkel'stein RG: 


Giamarchi-Schulz RG: 


disorder (perturbative) 


disorder (perturbative) 


& interaction (exact) 


& interaction (exact) 


Lt ~ C : strong 


Lt ~ ^ : strong coupling 


Altshuler-Aronov corrections 


in Giamarchi-Schulz RG 


7^ Anderson localization 


7^ Anderson localization 


~ 5, Affwi ~ CTD : 


~ ^, AcTwi ~ CTD : 


strong Anderson localization 


strong Anderson localization 


weak localization correction 


this work, 


to the conductivity, Eq. (|2.1[) 


Sees. IVIIVIII 


dephasing time, Eq. (|2.4|l 


this work. Sees. IVIIVIII 



order, described by the correlation function 

{U{x)U{x')) = S{x - x')/2tivqTq , (3.1) 

where tq is the elastic scattering time and vq — l/irvp is 
the density of states per spin in the absence of interaction. 
The Hamiltonian for spinful electrons takes the form: 



H — Hy_iu + Hcc + Hdis , 



(3.2) 



Hkin = '"Pif^k ~ kF)ipl^{k)ipi,^{k) , (3.3) 

k,fl — :L.(T 



Ha 



+ gl^lf,,^,^^,,a'^ , (3.4) 

(3.5) 



Here 



r{x) = ll^\,„{x)%k^,.a{x) 



(3.6) 



is the fermionic density, lAh and lAf are backscattering 
and forward-scattering amplitudes due to disorder. It 
can be shown22ii2i that the forward-scattering amplitude 
can be gauged out in the calculation of the conductivity. 
Therefore, from now on we set Uf = 0. The correlation 
functions for Ub{x) read 



{Ub{x)Ub{x')} 
{Ubix)U;{x')) 



, 

{U{x)U{x')) 



(3.7) 
(3.8) 



Two points are worth mentioning here. First, all the 
parameters of the above Hamiltonian include, in gen- 
eral, high-energy renormalization effects (similar to the 
Fermi-liquid renormalization in higher dimensionalities) 
and should be considered as effective couplings of the 
low-energy theory. Second, the low-energy theory is only 
then well-defined when supplemented by an ultraviolet 
energy cutoff A. The latter is determined by the Fermi 
energy, the bandwidth, or the spatial range of interac- 
tion in the original microscopic theory, whichever gives 
the smallest cutoff. 

Let us now concentrate on the spinless model. In 
this case, two essential simplifications occur. First, 
the interaction-induced backward scattering and forward 
scattering relate to each other as direct and exchange 
processes. Therefore, the backscattering Vh only appears 
in the combination V{ — Vh and thus merely redefines 
the parameters of the Luttinger model, formulated^i^ 
in terms of forward-scattering amplitudes only. In other 
words, the term with gi in Eq. (j3.4|) has the same struc- 
ture as the f;2-term, so that we can absorb gi into 52, 

52 = V{ - Vb = 32 - .91, 

and start our consideration with the following form of 
the Hamiltonian 



H^in = J2 ^F{fik-kF)4'l{k)i;^,{k) , (3.9) 
-^co = ^ ^ I dx{nf,g2n^f, + nf,g4nf,} ,(3.10) 
-ffdis = / dx{u^{x) 4>U- +Ub{x) (-3.11) 



In what follows we omit the tilde over the shifted cou- 
pling 92- The second simplification is related to the 34- 
term in Eq. (|3.10p . Naively, one could think it is simply 
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zero, since it describes a local interaction between identi- 
cal fermions. In fact, the situation is slightly more tricky. 
While yielding no genuine interaction, this term gener- 
ates a shift of the velocity due to a quantum anomaly, 

VF 



Vp, where 



Vp = Vp 



94_ 

27r 



(3.12) 



This is also straightforwardly seen in the bosonization 
approach, see Sec. ITVl 

Therefore, the complete set of parameters defining the 
spinless problem includes Vp, g2, tq, and A, with the 
interaction characterized by a single coupling g2. It is 
customary to define the corresponding dimensionless pa- 
rameter in the form 



K = 



1 - g2/2TTV*p 
1 + g2/2TTV*p 



1/2 



(3.13) 



Note that in a Id system with g2 ~ gA Eq- (13.13P is 
identical to another conventional representation of K in 
terms of the bare velocity vp, 



K = 



1 



(l+52/7r«f)i/2 



(3.14) 



For the case of a weak interaction (considered in the 
present paper), K ~ 1 — a, where 



52 



92 



2ttv*p 2-kvf 



< 1 



(3.15) 



is the dimensionless strength of the interaction. 



B. Qualitative discussion: Luttinger-liquid 
renormalization vs Fermi-liquid dephasing 



Having formulated the model in Sec. IIII Al wc turn 
now to a discussion of our strategy in its analysis. In 
principle, one could proceed in a number of alternative 
ways: 

(i) First of all, one can bosonize the action and try to 
study the transport properties of interacting electrons in 
the bosonic language. Sec. [iVl While this approach is 
very powerful for the problem of Id interacting fermions, 
it turns out to be poorly suited for the analysis of trans- 
port problems in the presence of disorder. For this rea- 
son, we develop an alternative "two-step procedure" . 

(ii) The two-step procedure combines the language of 
bosonization for high energies with a fermionic treatment 
of processes with low-energy transfers. This method is 
discussed on the qualitative level in the remainder of this 
subsection and is implemented in Sees. IIVIVII below. 

(iii) Finally, one more alternative is the "functional 
bosonization" approach This method also 
combines the languages of fermionic (electrons) and 
bosonic (plasmons) excitations, but in contrast to the 
two-step procedure this is done in a uniform fashion, so 



that the processes with high and low energy transfers are 
treated on the same footing. This method yields results 
identical to the two-step approach and is presented in 
SecEHl 

Let us describe now the ideas underlying our two-step 
approach. The method is based on the separation of two 
types of processes: (i) "elastic renormalization" of dis- 
order, which is associated with virtual transitions with 
energy transfers |w| > T and (ii) inelastic processes with 
< T which lead to dephasing in quantum-interference 
effects and thus to a destruction of the Anderson local- 
ization. This separation of energy scales is very much 
similar to that in physics of mesoscopic phenomena in 
higher-dimensional case, see Sec. Ill Al 



1. Step 1: Luttinger-liquid renormalization 

We begin by considering the Drude conductivity under 
the condition that the dephasing time is much shorter 
than the transport time of elastic scattering off disorder 
T. The source of the strong dephasing may be external 
(say, interaction with phonons, or Coulomb interaction 
with "environment") or intrinsic, that is the inelastic e-e 
interaction within the quantum wire. In the latter case, 
the limit <C r is achieved at sufficiently high T, as 
shown below in Sec. |Vll For <C r, the Anderson- 
localization effects are strongly suppressed by dephasing. 

To leading order in r^/r ^ 1, the conductivity is given 
by the Drude formula 



2 2 

e pVpT 



(3.16) 



where p is the compressibility. The Drude conductivity 
(|3.16p depends on T through a T-dependent renormal- 
ization of the static disorder i^iilliSiilSi^ 

t{T) = To(T/A)2- . (3.17) 

For spinless electrons, the exponent in Eq. (|3.17p reads 



a=l-K ' 



(3.18) 



In principle, the interplay of disorder and interaction 
leads to the renormalization of both disorder, Eq. (|3.17p , 
and interactionjSi However, the value of a in Eq. p.lSp 
is given by the bare interaction constant (the one in a 
clean system), see Sec. |TV]for details. 

The Luttinger-liquid renormalization l|3.17p is sim- 
ilar to the Altshuler-Aronov corrections^ [Eq. (|2.2p ] 
in higher dimensionalities. The underlying physics of 
the elastic renormalization of disorder can be described 
in terms of the T-dependent screening of individual 
impurities and scattering by slowly decaying in real 
space Friedel oscillations, similarly to higher-dimensional 
systemsi ^^i^^ At this level, the only peculiarity of Lut- 
tinger liquid as compared to higher dimensionalities is 
that the renormalization of r is more singular and ne- 
cessitates going beyond the Hartree-Fock (HE) approach 
even for the weak interaction/2i 
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Technically, Eq. (|3.17|) can be derived by a variety of 
methods: either fermionic^i^^i or bosonic j?^'"? In partic- 
ular, a conceptually important framework is based on 
the bosonic RG procedure developed by Giamarchi and 
SchulzJi This approach allows one to integrate out all 
energy scales between A and T, which are responsible for 
the elastic Luttinger-liquid renormalization due to virtual 
e-e scattering processes. We will discuss the features of 
the bosonic RG approach^ relevant to our analysis in 
more detail in Sec. HVl 

The renormalization of r stops with decreasing T at 
Tt(T) ~ 1. Physically, this is due to the fact that the 
long-range Friedel oscillations are cut off even at zero T 
on the spatial scale of the order of the disorder-induced 
mean free path. The condition 

Tt{T) - 1 (3.19) 

gives the zero-T localization length 

e oc To^/(i+2'^) . (3.20) 

It is important to stress, however, that Eq. (|3.19p does 
not correctly predict the onset of localization, in con- 
trast to the argument made in Refs. [33l8ll . This can 
be seen, in particular, by noting that the temperature 
T ~ does not depend on the strength of interaction 
for small a, whereas it is evident that for noninteract- 
ing electrons there is no interaction-induced dephasing 
(the dephasing length is infinite, = oo) and hence 
cr(T) = for any T due to the Anderson localization, 
see Sec. Ill B II The onset of strong localization occurs 
at ^ ^, as in the case of higher dimensionality, see 
the discussion in Sec. Ill Al The temperature Ti below 
which the Anderson-localization effects become strong is 
determined in Id by the condition 

T(ri) ~ T^(ri) . (3.21) 

At this temperature, the quantum-interference effects 
leading to the WL correction ActwI to the conductivity 
are no longer weak, 

|Aawi(Ti)| ^t7D(Ti) . (3.22) 

2. Step 2: Fermi-liquid dephasing 

It is worth stressing that the Luttinger-liquid physics in 
the present problem is fully accounted for by Eq. (I3.17p . 
i.e., all the Luttinger-liquid power-law singularities are 
now incorporated in the renormalization of t(T). At this 
point, we are left with the system of interacting electrons 
with a new effective "Fermi energy" T. These electrons 
are scattered by disorder with renormalized tq — > t{T), 
which is still weak for Tt{T) > 1. 

What remains after the renormalization is the e-e scat- 
tering with energy transfers \uj\ ^ T. The interaction is 
thus still important and leads to real inelastic processes. 



The physics of these processes is in essence equivalent to 
the Fermi-liquid physics of dephasing in higher dimen- 
sionalities. Having integrated out the virtual transition, 
we have to choose a controllable method of dealing with 
the dephasing. The bosonization approach (whose ad- 
vantage is that it can be straightforwardly formulated 
for the arbitrary strength of interaction) is not particu- 
larly beneficial here. Although, in principle, it contains 
all effects of the Anderson localization in Id, the ma- 
chinery of bosonization is poorly suited to deal with such 
ingredients of the localization theory as the interference 
and dephasing. Even for noninteracting electrons, it is 
by far not straightforward to derive the Anderson lo- 
calization [for instance, to obtain the Berezinskii-Mott 
formula (|2.6p for the case of weak disorder] in the lan- 
guage of bosonization. For this reason, we refermionize 
the theory at this step. Sees. IVIVII 



IV. BOSONIZATION APPROACH 
A. Bosonized action 

As discussed in Sec. lIIIBl the first part of our program, 
namely the renormalization of the model due to high- 
energy processes, is done most efficiently in the frame- 
work of the bosonization technique. In this section, we 
present a calculation of t{T), using the RG scheme devel- 
oped by Giamarchi and Schulzi^i After a brief description 
of the derivation of the RG equations, we will discuss an 
important subtlety of the scheme,™ namely a mixing of 
disorder and interaction. 

For simplicity, we consider here the spinless case (in 
notation of Ref. we set gi± = g2± = and hence 
Up = Uc, Kp — K„). As discussed in Sec. IIII Al for 
spinless electrons the backscattering amplitude gi can 
be absorbed into 172, and 174 into v*p, so that only one 
interaction coupling constant remains (32), yielding the 
dimensionless Luttinger-liquid parameter K, Eq. (|3.13[) . 

Using the boson representation of the fermionic 
operators, the Hamiltonian can be expressed in terms 
of the boson operators (/> and 11 that satisfy the com- 
mutation relation [n(x), 4>{x')\ — —i5(x — x'). In a given 
realization of disorder the Hamiltonian reads 

H = dx[uK{'Knf + ^{d,^f] 

+ ^y"da;{Wb*(x)exp(2i0) -f H.c.} (4.1) 

with 

u=[{vlf^gl/A.'f' . (4.2) 

The ultraviolet length scale A is related to the energy 
cutoff A as 

A = u/nK . (4.3) 
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We will concentrate on the case of e-e repulsion, implying 
that K < 1. The disorder-induced forward scattering 
amplitude Uf{x) is absorbed^ in the quadratic part of 
H as a shift of <j>{x). 

In order to perform the disorder averaging, one intro- 
duces replicas, 0„. In the resulting action the interac- 
tion (Luttinger- liquid) term Sll is quadratic in fields 0„, 
while disorder generates a term Sd of cos 20-type which 
breaks the Gaussian character of the theory: 



(4.4) 



Sll 



dxdT 
uK J 27r 



(4.5) 



™ dx c?T cZt' 

-Db / ^TT-TT^ cos{2[0„(a;,r) - 0™(a;,T')]} 



(27rA)2 



(4.6) 



where Di, — Vp/2tq. The integration over the imaginary 
times is performed within the interval < r < 1/T. 

The conductivity is expressed in terms of the current- 
current correlation function at zero momentum via the 
Kubo formula: 



with 



e 

in 



I" Jo 



/C(r) - / dx{Trj{x,T)j{0,0)) , 



(4.7) 



(4.8) 



where in the bosonic language the currents are expressed 
via the field (t>(x, r) as 



j{x,T) = -dr(f'ix,T) . 



(4.9) 



The first (diamagnetic) term in Eq. (|4.7p is canceled by 
a contribution of the singularity at r = in the sec- 
ond term. It is worth noting that, strictly speaking, the 
correlation function {Trdr4>{x,T)dr4>{(^T^)) is a reducible 
(with respect to the interaction) correlator, when it is 
calculated with the use of Sll as given by Eq. (|4.5p . On 
the contrary, the conductivity is given by the irreducible 
current-current correlator. However, at zero external mo- 
mentum Q, the two correlators are equal to each other, 
since the polarization operator vanishes at Q = for 
arbitrary nonzero frequency. This should be contrasted 
with a calculation of the conductance of a finite inter- 
acting wire: the latter is related to the response at finite 
Q and the Kubo formula should be supplemented with 
an analysis of the effect of renormalization of the electric 
field by e-e interaction i^^'^^ 

Note that the authors of Ref. [13 did not calculate the 
conductivity using the Kubo formula (|i?7)) . 
Instead, they found the renormalized value of l/r and 
substituted it into the Drude formula, Eq. p.l6p (which 



can be represented as a certain approximation in the 
framework of the memory- function formalism, see, e.g., 
Refs. [78ll97[ ). In view of technical difficulties, in this paper 
we do not attempt to evaluate the conductivity, including 
the localization effects, by using the Kubo formula in the 
bosonic language either. Having used the bosonization 
for the analysis of the high-energy renormalizations, we 
refermionize the theory and calculate the conductivity in 
the fermionic language. Sec. IVII 



B. Giamarchi-Schulz RG 

Giamarchi and Schula^l performed the renormalization 
of the coupling constants u, and Di, in the action (14. 4p 
upon rescaling the ultraviolet cutoff X ^ L = Aexp(£). 
Their analysis is done to first order in the dimensionless 
strength of disorder 



(4.10) 



To this order, one can omit replica indices in Eq. (|4.4p . 
Before presenting the RG equations, we discuss a subtle 
point of the approach developed in Ref. [H. Specifically, 
the coupling constant K in terms of which the RG equa- 
tions are obtained in a natural way, is not simply K but 
rather contains an admixture of disorder. We now briefly 
reproduce the corresponding argument from Ref. [s^ . In 
the fermionic representation, the part of the action re- 
lated to the disorder- induced backscattering acquires, af- 
ter the disorder averaging, the form 



Db I dx dr dr' tjj^^{x,T)ip-^{x,T) 



Sd — I Lt.cu,/ u,/ yj^ 

X ijUx,T')'ilj+{x,T') 



(4.11) 



In order to derive the RG equations, Giamarchi and 
Schulz singled out the ultraviolet contribution of close 
times u\t — t'\ < A, splitting the time integration into 
two parts: 

Sd — ^Db / dxdrdr' 

J u\r — r'\>\ 

X ^l{x,T)^.{x,T)i^l{xy)^+{x,T') (4.12) 
dx dr 



(4.13) 



X ■4;\_{x,T)i;_{x,T)'4>'^_{x,T)'4>+{x,T) 



The second term in Eq. ()4.13|) is local in time and 
thus equivalent to an interaction term generated by 51- 
processes. For this reason, the authors of Ref. rede- 
fined the interaction constant, gi = gi — 2DbX/u. For the 
case of spinless electrons considered here, this is equiva- 
lent to a modification of the constant g2- 



92 ^92 + 2DbX/u 



(4.14) 



10 



In terms of the modified interaction constants the RG 
equations^ (in the spinless case only three equations out 
of original six remain) are: 



dk/di = , 

du/M = -uKV/2 , 
dV/di = {S-2K)V 



(4.15) 
(4.16) 
(4.17) 



where u and K are related to the modified interaction 
constant 



(4.18) 



according to Eqs. (|4.2p and (|3.13p . respectively. 

The modification of the plasmon velocity u by dis- 
order, as well as its renormalization, can be neglected 
in Eq. (|4.18p . since this generates higher powers of V. 
Then Eq. (|4.16p . which describes the renormalization of 
M, turns out to have no influence on the renormalization 
of disorder and interaction, described by the remaining 
two equations. [In fact, Eqs. (|4.15p . (|4.16p show that K/u 
is the invariant of the flow.] Therefore, we will concen- 
trate on the coupled set of the flow equations (|4.15p . 
(I4T7)) . It follows from Eq. (|4T7)) that for K < 3/2 
the strength of disorder ^{i) increases upon renormal- 
ization. At the same time, the value of K decreases. The 
system flows then to the strong-coupling regime where 
P becomes of order one and the first-order treatment of 
disorder is no longer valid. 

The decrease of K{£) upon renormalization suggests 
that the initially weak interaction {K close to 1) becomes 
strong upon approaching the strong coupling regime. 
However, this is not so, as we demonstrate below: the 
decrease of K occurs solely due to the admixture of in- 
creasing disorder in this quantity. 

It is also worth noticing that Eq. (|4.15p seems to gener- 
ate finite interaction due to the scattering off disorder in 
the noninteracting case. However, as verified in Ref. [s^, 
this is not the case. The point is that the condition K = 1 
does not mean that the system is noninteracting, because 
of the mixing of interaction and disorder in Eq. (|4.18p . 
Therefore, it is desirable to rewrite Eqs. (|4.15p . (|4.17p 
in terms of the true interaction constant K instead of 
the one modified by the admixture of disorder, K (see 
Fig. 9.6 in Ref.©. 

Let us first consider the case of weak interaction, 1 — 
ii' ~ a <c 1, retaining only the terms of first order in a 
in the RG equations. From Eqs. (jXT^ . fO)) . (liHI)) . and 
(I4.14P we get, to first order in a, 



k = l-a-V/2 + aD/2 



(4.19) 



Substituting this formula in Eq. (|4.17p . we express the 
flow equation for the disorder strength in terms of the 
renormalized true interaction constant a{t): 



In this equation, the first term in the right-hand side 
describes Ohm's law, while the second term (oc a) yields 
the Luttinger-liquid renormalization of r, Eq. (|3.17p . 

The last (quadratic-in-disorder) term in Eq. (|4.20p 
is beyond the accuracy of the derivation of the RG 
equations and should be neglected. We note, how- 
ever, that thi s sp urious term is hidden in the RG equa- 
tions of Ref. [SJ. Accidentally, this spurious quadratic 
term qualitatively mimics the effect typical for the WL. 
Clearly, this purely accidental similarity does not hold 
on the quantitative level: the true WL — when included 
into consideration — would only appear in the r.h.s. of 
Eq. lOO)) at X»3 order. 

Now let us turn to the flow equation (|4.15p for the 
interaction strength. We rewrite it using Eqs. (I4.19P and 
in the form 



da/dl = -3aX'/2 - 



(4.21) 



Again, the last (quadratic-in-disorder) term here is be- 
yond the accuracy of the calculation and has to be ne- 
glected (this unphysical term generates interaction from 
disorder) . It is important that the linear-in-disorder term 
in Eq. (|4.2ip vanishes for a = 0. In other words, disor- 
der does not generate interaction, which was not trivial 
to see from Eq. (|4.15p . 

To summarize, for a weak interaction, a <C 1, the set 
of two scaling equations for a and T) reads 



da/dl = -3aX'/2 , 
dV/d£ = (1 -I- 2a)V . 



(4.22) 
(4.23) 



One sees that Eq. (|4.22p yields decreasing a{£): while 
the disorder strength grows with increasing L, the inter- 
action gets weaker ("disorder kills interaction" — Ref. HI; 
also Ref. l98h . in contrast to a naive expectation based 
on Eq. (|4.15p . This resembles the behavior of the renor- 
malized interaction in 2d diffusive systems. Specifically, 
according to the Finkel'stcin RG equations, the combi- 
nation of the singlet and triplet amplitudes (an analog of 
a for the spinful case) which determines the interaction- 
induced correction to the conductivity Actaa, decreases 
for weak interaction (a <C 1) due to the disorder- induced 
renormalization of interaction. 

For completeness, let us also consider the RG equations 
for arbitrary strength of interaction. Expanding K in'D 
and using Eqs. (|4.2p and (|4.14p . we obtain 



k ~ K - {1 + K^)V/A 



(4.24) 



Substituting Eq. in Eqs. ((i7T5)) and (|iTfp . we ar- 

rive at 



dK 

Hi 



V , (4.25) 



(4.26) 



dV/de^V + 2aV + V^ 



(4.20) 



The solution to Eqs. (|4.25p . (|4.26p is shown in 
Fig. [1] for comparison, we show also the solution to 
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0.25 0.5 0.75 1 1.25 1.5 1.75 



K 

FIG. 1: Upper panel: the RG flows given by 
Eqs. (|4.25p . (|4.26|l in the plane {K,!)), where K is the true 
interaction constant and T> is the strength of disorder (the di- 
mensionless resistance). Dashed line: the separatrix touches 
the quantum critical point a.t K = 3/2, 2? = 0. Long-dashed 
line: the line of if = 1 separates the repulsion [K < 1) and 
attraction (K > 1) regions. Lower panel: the RG flows from 
Eqs. (|4.25p . (|4.26|l in the coordinates K and V/L, the latter 
corresponds to the resistivity. The renormalization for each 
curve is stopped at D = 1 (the strong coupling regime). The 
dashed line separates the repulsion and attraction regions: 
an attractive (repulsive) interaction cannot become repulsive 
(attractive) in the course of the renormalization. Initial cou- 
plings in the upper panel from left to right: Do = 0.001, 
Ko=0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3 
and Ko = 1.9, ©o = 0.3, 0.25, 0.2, 0.15, 0.05. Sep- 
aratrix: Ko ~ 1.9; Do ~ 0.079. Initial couplings in 
the lower panel: A" < 1 (from left to right): Vo — 
0.001, Ko = 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9; 
K > 1 (from top to bottom): Ko ~ 1.9, Vo — 
0.4, 0.35, 0.3, 0.25, 0.2, 0.15, 0.10, 0.05. 



Eqs. (liTSjl . fCTT)) in Fig. H The flow equation 
for K tells us that the repulsive interaction, K < 1, 
gets weaker in the course of renormalization, dK/d£ > 0, 
whatever the interaction strength is, see Fig. [TJ Equa- 
tions (|125l),(|12n]) reduce to the set of Eqs. gH]) , (gini) 
around K = 3/2, where the Giamarchi-Schulz RG pre- 
dicts the metal- insula! or transition. Note that for K > 
3/2 the flows in terms of K and K are qualitatively sim- 
ilar to each other, as is evident from a comparison of 
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0.25 0.5 0.75 1 1.25 1.5 1.75 



K 

FIG. 2: Upper panel: the RG flows given by 
Eqs. (|4.15p . (|4.17|l in the plane of the modified interaction con- 
stant K, introduced in Ref. [s^ . and the strength of disorder 
(the dimensionless resistance) D. Dashed line: the separa- 
trix touches the quantum critical point nt K = 3/2, V = 
0. Lower panel: the RG flows from Eqs. (|4.15p . (|4.17|l in 
the coordinates K and T>/L (corresponding to the resis- 
tivity). The renormalization for each curve is stopped at 
T> = 1 (the strong coupling regime). Initial couplings in 
the upper panel from left to right: = 0.001, Ko = 

0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3 
and ko = 1.9, Vo = 0.3,0.25,0.2,0.15,0.05. Separatrix: 
Ko = 1.9; Vo — 0.104. Initial couplings in the lower panel: 
Vo = 0.001, ko = 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 and 
ko = 1.9, Vo = 0.4, 0.35, 0.3, 0.25, 0.2, 0.15, 0.10, 0.05. 

Figs. □ and H 

Let us return to the case of weak interaction and an- 
alyze the RG flow generated by Eqs. ((1221) , (|1!2S1) , see 
Fig. [3l An approximate solution of these equations in 
the range of weak disorder, V^f) <C 1 (which is where 
they are valid), is easily found by an iterative procedure. 
The result is: 

I?(^) ~ Po , (4.27) 

a(£)~ao[l-3I?(^)/2] , (4.28) 

where I?o and ao are the initial (ultraviolet) values of the 
coupling constants at L ~ A. 

It is worth emphasizing that whereas the correction to 
q;o is of the order of ao itself when 'D{(.) ~ 1, the renor- 
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FIG. 3: Upper panel: the renormalization of the interac- 
tion constant K due to disorder for various initial values of 
K in the ultraviolet limit A/L = 1. The strength of interac- 
tion decreases upon renormalization, K ~* 1. Lower panel; 
the renormalization of the "conductivity" L/T) in units of the 
bare "conductivity" A/Do (solid line). The dashed lines de- 
scribe the dependence of L/T) given by Eq. (|4.27|l with the 
exponent ao = 1 — Kq, where Kq is the ultraviolet value of 
the Luttinger constant. Each flow is stopped at D = 1. Ini- 
tial couplings from top to bottom, both panels: 23o = 0.001, 
Ko = 0.95, 0.9, 0.85, 0.8, 0.75, 0.7, 0.65, 0.6, 0.55, 0.5. 



malization of disorder is governed by the initial value of 
the interaction constant and not by the running coupling 
constant a{C). The renormalization of the interaction 
constant is irrelevant in the present case, and hence the 
renormalization of disorder reduces to the renormaliza- 
tion of an individual impurity. As a consequence, the 
exponent a in Eq. (I3.17|) for the renormalized scatter- 
ing time t{T) is given by the bare interaction constant 
(the one in a clean system), a = ag ~ ao. Thus, in 
the Drude formula (|3.16p the exponent of the Luttinger- 
liquid power-law function (T/A)° is T-indcpendent. 



C. Bosonization vs Localization 

As discussed in Sec. lIVI3| the solution of the RG equa- 
tions flows to the strong-coupling regime, where V > 1 
and the perturbative-in-disorder RG equations are no 



longer valid. Since in the absence of interaction the 
strong-coupling fixed point I?/L — > oo corresponds to the 
complete Anderson localization, Giamarchi and Schulz 
identified the regime I? ^ 1 as the localized phase. 

It should be stressed, however, that the notion of An- 
derson localization does not follow from the RG equations 
of Ref. fsil. Indeed, Eq. (|4.23p does not account for the 
renormalization of disorder V/ L in the absence of interac- 
tion: the effects of Anderson localization (arising due to 
interference terms which involve higher powers of V) are 
not included in this RG scheme. In the weak-coupling 
regime, I? <C 1, the Anderson localization would show 
up in an additional a-independent term, fw\{'D) ~ 
on the right-hand side of Eq. (|4.23[) , which is beyond the 
accuracy of the derivation of these equations. It is also 
worth noting that in the presence of interaction all the 
RG trajectories for V/L tend to saturate at finite val- 
ues of V/L, as can be seen from Fig. [1] in contrast to 
what one expects for the Anderson-localization behav- 
ior. The strong-coupling fixed point I? ^ oo is therefore 
a manifestation of Ohm's law, rather than of the Ander- 
son localization. 

In order to study the temperature behavior of cr(T), 
the renormalization procedure in Ref. was stopped at 
L = Lt ^ ti/T, with the argument that the coherence 
effects disappear on length scales greater than Lt due 
to inelastic processes. Although this is perfectly correct 
for the interaction-induced renormalization of disorder 
which is cut off by the thermal smearing of the electron 
distribution function [similarly to the Altshuler-Aronov 
correction to the conductivity in higher-d systems, see 
Eq. (|2.2p ]. the coherence effects relevant to the Anderson 
localization are governed by the dephasing length, L^. 
The condition Lt ~ C yielding the temperature ~ u/^ 
for the onset of the Anderson localization is in general 
not correct. 

Indeed, the additional term f^\{T>) responsible for the 
WL (quantum-interference) effects is cut off not by Lt 
but rather by the dephasing length L^, which is paramet- 
rically different from Lt- In particular, in the absence 
of interaction L^ diverges, while Lt remains finite, be- 
ing independent of a. For weak interaction L^p 3> Lt 
and hence, there exists an intermediate range of length 
scales, Lt ^ L < L^, in which only the WL term fw\{V) 
contributes to the renormalization of V. Thus, it turns 
out that the scalings with T and L are not interchange- 
able, when a{T) is considered. For a ^ 1, the onset of 
strong localization is determined by the condition 

i0 ~ e , (4.29) 

which in Id is equivalent to 

~ T . (4.30) 

After the RG transformation, the bosonization action 
retains its original form (|4.4p - (|4.6p . but with a renormal- 
ized ultraviolet cutoff, A ^ T, and a renormalized scat- 
tering rate, Tq^ — > t~^{T). Performing a transformation 
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from the bosons back to fermions, we thus return to the 
original fermionic theory with a corresponding renormal- 
ization of its parameters. This theory is analyzed in the 
following sections. 



V. GOLDEN RULE IN LUTTINGER LIQUID 

In order to analyze the effects of dephasing, we will 
follow the route suggested by earlier works on higher- 
dimensional (in particular, diffusive) systems. Since the 
physics of dephasing is governed by electronic inelas- 
tic scattering processes, a natural first step is to cal- 
culate the e-e scattering rate t^J^. Indeed, for higher- 
dimensional systems in the high-temperature (ballistic) 
regime, Tt ^ 1, the dephasing rate t^^ to leading order 

is given by t~^, see Ref-lool For quasi-ld and 2d systems 
in the diffusive regime the scattering kernel acquires an 
infrared singularity leading to a divergent t~^. How- 
ever, even in this case the calculation of the e-e collision 
rate turns out to be instructive: a parametrically correct 
result for the dephasing rate can be obtained from the 
expression for t^^^ supplemented with an appropriate in- 
frared cutoff. It is thus useful to begin by analyzing t^^^. 
We will show that this quantity is finite and meaningful 
in the interacting Id system (contrary to what one might 
expect based on the famous non- Fermi-liquid character of 
the Luttinger liquid). On the other hand, we will demon- 
strate that, while governing the dephasing of Aharonov- 
Bohm oscillations, does not give the dephasing rate 
for the WL of spin-polarized electrons. 

The Golden Rule expression for the e-e collision rate 
following from the Boltzmann kinetic equation reads 



diij J de' IC{uj) 



Too(e) 

X {ftj.'fe'+^+h-uf^.h'+u) , (5.1) 

where /C(w) is the kernel of the e-e collision integral 
St{/J = jdio j de'JC{oj) 

X (~/c/c-w/c'/e'+w ft fe—i^fe'ff-'+i^) 



+ I doj I de'IC{uj)f,.^f^,f,, 



Tcc(e) 



+U1 



Here is the Fermi distribution function, f^ = l — f^. 

To leading (second) order in interaction, the scattering 
kernel K{uj) in a clean Luttinger liquid has the form 

K{uj) = 7^,[Ki^+{uj) + K^_{uj)] + K^{lu) , (5.2) 

where 

= ■94^/J[R^^+(-''^)]' ' (5-3) 
= gl f p-ReD+{Lu,q)ReD-iuj,q) , 



= -K 



H 



(5.4) 
(5.5) 



Here K^j^ and are related to scattering of two elec- 
trons from the same or different chiral spectral branches, 
respectively, Kp is the exchange counterpart of 
r]s is the spin degeneracy (77s = 1 for the spinless case 
and rjs — 2 for the spinful case), and D±{uj,q) are 
the two-particle propagators, given by Eq. (jA.4p with 

Substituting Eqs. ((0)) - (|53)) in Eq. ((ET|) . we obtain 
the lowest-order result for the e-e scattering rate at the 
Fermi level (e = 0) in terms of the corresponding contri- 
butions to the retarded electronic self-energy E;^ defined 
by Gl{e,p) = [e - vpp - I]^(e,p)]~\ where is the 
retarded Green's function for right-movers: 



l/2ree = -ImE^(0,0) , 



(5.6) 
(5.7) 



The diagrams for and are shown in 

Figs. UK, HId, and [H;, respectively. The contribution of 
the diagram in Fig. [Ha (Hartree term) is canceled by the 
diagram in Fig.[4j: (Fock term) due to the Pauli principle. 



(5.8) 



We refer to this cancellation as the "Hartree-Fock can- 
cellation" . 



a) 






FIG. 4: Second-order ("Golden Rule") diagrams for the self- 
energy of a right-mover. Wavy lines: interaction of the right- 
mover with left- ((72) or right- (34) movers, (a) the 
Hartree contribution due to the interaction with a left mover; 
(b) , the Hartree contribution due to the interaction with 
a right mover; (c) E^, the Fock (exchange) contribution. 



Let us now calculate the imaginary parts of the per- 
turbative self-energies. Using Eq. (|A.4p to get the re- 
tarded propagators D±(uj,q) by analytically continuing 
iQn —^uj + iO and putting 52 = .94 = 27rwi?Q; we find: 



ImE 



H 



du ( 
2-K \ 



coth tanh — ) 

2T 2T) 



ZTT 



diij I dqS{uj — VFq)S{iLj ~ vpq) 



(5.9) 
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jT ^ ^ „ f duj / 10 u) \ 

ImS^ =-2t\^olv%\ — w coth tanh — 

X [ — RcD+{Lj,q) ReL>_(w,g) ~ -na^vpT 
J 27r 



duj 

\<T 



dqS{uj — vpq)S{uj + vpq) 



(5.10) 



One sees that the contribution of is diverging. For 
spin-polarized electrons (i.e., for 77.; = 1) it is, however, 
exactly canceled by the exchange interaction, Eq. (|5.8p . 
Indeed, as we have discussed in Sec. IIII Al the 34 interac- 
tion drops out of the problem in this case, inducing only 
a shift of the velocity, Eq. (|3.12p . The remaining term 

-!■ and is given by 



S:^_ is determined by us, q 



ImE 



H 



(5.11) 



Already the perturbative expansion demonstrates a qual- 
itative difference between two cases of spinless and spinful 
electrons — the spinful case will be considered in Ref. 0. 

The characteristic frequencies in Eqs. (|5.9p . (|5.10p sat- 
isfy \uj\ < T, because of the factor coth(a;/2T) — 
tanh(a'/2T), even if we would perform this calculation 
in the original theory with the ultraviolet cutoff A. For 
the sake of comparison, we present also the corresponding 
contribution to the real part of the self-energy: 



X < coth — 



^ -Siuj + ypq) 



tanh 



2T e + LO - vf{p + q) 
uj + e q 



2T uj + ypq 



S[e + Lu -vpip + q) _ 



which for T ^ \e ± vfp\ gives 



IteE+((,p) 



-ion. 



VFP) In ^ 



(5.12) 



In contrast to the imaginary part of the self-energy, the 
real part is determined by energy transfers \lu\ ^ T, which 
is characteristic of the elastic virtual processes. Further- 
more, Rel]+ is logarithmically divergent in the ultraviolet 
and is cut off by A: this behavior is specific to Id. 

It is worth stressing that in the considered case of weak 
interaction, a <C 1, the lowest-order contribution to the 
inelastic scattering rate obtained from Eqs. (I5.6p . (|5.1ip . 



na^T 



is much smaller than temperature: 

Tree > 1 . 



(5.13) 



(5.14) 



In higher dimensionalities, this fact is commonly referred 
to as one of the conditions of the existence of Fermi liq- 
uid: the inverse lifetime of a quasiparticle is smaller than 
its characteristic energy. In this respect, the weakly in- 
teracting Luttinger liquid, while being a canonical exam- 
ple of a non-Fermi liquid, reveals the typical Fermi-liquid 



property. The Luttinger-liquid physics is in fact encoded 
in the singular real part of the perturbative self-energy, 
Eq. ([5T^ . For a^ln{A/T) < 1, the Green's function 
can be written as G^{e,p) ~ Z/ {e — vpp + ma^T/2), 
where the single-particle weight Z ~ 1 — [0^/2) ln(A/T) 
is suppressed by interaction. It is in fact the product 
ln(A/r), coming from the real part of the self-energy, 
that after resummation of all orders gets exponentiated 
and leads to the Luttinger-liquid power-law singularities. 
Those singularities have been accounted for by the renor- 
malization at step 1 in the present approach. 

It is instructive to compare the perturbative collision 
rate in a clean Luttinger liquid, Eq. (|5.13p . with the 
dampingiSS of the exact retarded single-particle Green's 
hmction in the coordinate-time (x, t) representation. 



gl{x,t) ^2iQ{t)lin-g±(x,t) 



(5.15) 



Considering right-movers for definiteness, we have 

g+{x,t) = (A/27ru)Ai+"''/2(a:/u-t + i/A) 

X A^"'"^ {x/u + t-i/K) , (5.16) 

where 

A{z) = ttT/A sinh(7rrz) (5.17) 

and ah = [K + K^^ — 2)/2, which becomes ah — o? 12 
for a <C 1. 

The Green's function is peaked at the "classical 
trajectory" a; = ut. In the limit of large t ^ 00 we have 
an exponential decay^ 



, -iA /27rT\"''/' , 

g+{x ^ ut,t) ~ — — exp -TrafeTi) , 

2ttu \ lA J 

(5.18) 

which should be contrasted with a power-law t- 
dependence at zero temperature. 



g+{ut,t) = 



-iA 



1 



2tiu (2iAt)"''/2 



(5.19) 



The large factor A in the numerator of Eq. (|5.18p is remi- 
niscent of the pole in the noninteracting Green's function. 
It is worth stressing that, since |v4(i/A)| 1, the non- 
analytical singularity of the factor sinh~"''^^[7rT(x/w — 
t + i/A) ] in the main peak a,t x = ut is exactly canceled 
by the Luttinger-liquid rcnormalization factor (r/A)"''/^ 
which arises due to the resummation of the real part of 
the self-energy. Thus the exponential decay in Eq. ()5.18p 
comes solely from the factor [A{2t — z/A)]"''/^ which de- 
scribes the tail of the "left-moving peak" centered at 
X = —ut. The decay exp(— TrabTt) of the "residue" 

[x-u{t- i/A)]gf{x, t%^ut oc sinh-"''/2(2^Tt) 



for i ^ 00 agrees with the Golden Rule expression (|5.13p 
to order a^: irabT ~ l/2Teo. Therefore, the temporal de- 
cay of the Green's function is determined by the processes 
of inelastic scattering of right-movers on left-movers. 
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We see that the perturbative result for the scattering 
rate yields the asymptotic behavior of the exact Green's 
function. This means that the Golden Rule calculation 
of the inelastic scattering time is meaningful in Luttinger 
liquid. However, the Golden Rule time Tcc does not cut 
off the WL correction in the spinless case, as we demon- 
strate below in Sec. IVII 



VI. WEAK LOCALIZATION AND DEPHASING: 
FERMIONIC PATH-INTEGRAL APPROACH 

A. WL dephasing vs Golden Rule 

The notion of dephasing associated with the behav- 
ior of the single-particle Green's function (|5.20p makes 
sense in a clean Luttinger liquid in the ring geometry, 
see Sec. IVIEI However, as we are going to show in this 
section, the WL dephasing time r^' that governs the cut- 
off of the WL correction Auwi is parametrically different 
from the Golden Rule time Too, Eq- (|5.13p . 

This difference can in fact be anticipated based on the 
following observation. The integral (|5.10p for the Golden 
Rule rate is determined by the energy and momentum 
transfer q,uj = 0. On the other hand, in the spirit of 
Ref. [HI, soft inelastic scattering with qvp,^ ^ l/'''^' is 
not expected to contribute to the WL dephasing. This 
suggests that the dephasing rate l/r^' requires a self- 
consistent cutoff in Eq. (|5.ip at qvp,uj ~ ^/t^^, similarly 
to Eq. (|2.4p for higher-dimensional systems, and so is 
parametrically different from the one in Eq. (j5.13p . In 
particular, if one introduces, similarly to higher dimen- 
sionalities, a sharp low-frequency cutoff in the right-hand 
side of Eq. (|5.10p . one immediately gets zero instead of 
Eq. (|5.13p . On the other hand, introducing disorder in 
D+ and H_ in Eq. (|5.10p . one effectively broadens the 
delta functions by l/r, and the self-consistent equation 
can be solved. 

Such a self-consistent equation should be treated with 
the utmost caution, in view of the singular character 
of the ballistic propagators. In particular, when trans- 
formed to the space-time representation, the integral in 
Eq. (|5.10p is determined by a single point, x,t = 0, so 
that one could be drawn to the erroneous conclusion that 
it is high — rather than low — frequencies that matter. In 
any case, this equation is at most qualitative. However, 
the message that we learn from this consideration is cor- 
rect: (i) the WL dephasing rate differs from the Golden 
Rule result, and (ii) to evaluate it, we have to include 
disorder in the interaction propagators. To perform a 
controllable quantitative analysis, we employ the path- 
integral technique, which is a natural tool in the situa- 
tion when the diagrammatic calculation suffers from the 
infrared problems. 



B. Weak localization in Id 

We begin by identifying the leading contribution to 
the WL correction to the conductivity in the regime of 
strong dephasing, t'^^/t <C 1. This regime takes place 
at sufficiently high temperatures — the precise condition 
is found below in Sec. IIIc, where we evaluate tJ^. The 
diagrams that give the main contribution to ActwI are 
shown in Fig. [S] 




b) c) 

FIG. 5: Diagrams describing the leading WL correction to 
the conductivity of Luttinger liquid for r"' <g r. The dashed 
lines represent backscattering off impurities. The current ver- 
tices are dressed by impurity ladders. The solid lines are the 
Green's functions with disorder-induced self-energies. The di- 
agrams are understood as fully dressed by e-e interactions. 

The "three-impurity Cooperon" in Fig. (5^- describes 
the propagation of two electron waves along the path 
connecting three impurities ("minimal loop") in time- 
reversed directions. Fig. [6^. Indeed, one can easily check 
that the interference paths involving only two backscat- 
tering impurities are impossible merely due to geomet- 
rical reasons. This means that the Cooperon with two 
backscattering impurity lines^'^-'- yields zero contribution 
to the conductivity to leading order in (ei?r)~^. If one 
retains the forward scattering and does not gauge it out, 
the two-impurity Cooperon diagram (Fig. [7^) is canceled 
by the contribution of two other non-Drude diagrams 
with two impurity lines (Fig. [7t),c). The diagrams rep- 
resented in Fig. [5)3, c (the corresponding trajectories are 
shown in Fig. [61d) are of the same order as the Cooperon 
diagram. 

In the absence of interaction, all higher-order di- 
agrams (including non-Cooperon diagrams) which de- 
scribe quantum-interference processes involving a larger 
number of impurities (4, 5, . . .) are equally important and 
together with the diagrams in Fig. [5] sum up to ex- 
actly cancel ctd, which spells a complete localization.2^ 
In the WL regime, they only yield subleading corrections 
through a systematic expansion in powers of the small 
parameter tJ^/t. Thus, the WL in Id is governed by 
the "shortest possible Cooperon" , namely, the Cooperon 
with three impurity lines. In essence, this is due to the 
fact that the WL condition 1/t^' ^ reduces in Id to 
the condition of the "ballistic dephasing", tJ^/t ^ 1. In 
turn, this means that the contribution to Acr^i of typical 
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FIG. 6: (a) Interference paths described by the Cooperon 
diagram, Fig. [S^. Black dots: the backscattering off impuri- 
ties. The one-dimensional graphs are stretched vertically for 
ease of visualization. The two waves start at point i and in- 
terfere at point / after traversing the closed loop in opposite 
directions, (b) Interfering paths for the diagram in Fig. [5};. 
The dashed circle denotes the absence of backscattering off 
the impurity at point Xc- The paths corresponding to the 
diagram in Fig. [SJj are obtained by interchanging the upper 
and lower plots in Fig. |6Jd. Note that similar paths describe 
the WL in the ballistic regime in higher dimensionalities^^ 




FIG. 7; Two-impurity non-Drude diagrams. 



configurations of three impurities separated by a distance 
of order vpr is strongly suppressed by dephasing. In 
the compact configurations that do contribute the three 
impurities are anomalously close to each other, with a 
characteristic distance between them being much smaller 
than the mean free path. 



C. WL dephasing in Id: Qualitative discussion of 
the leading-order cancellation 

Before turning to the path-integral evaluation of the 
WL dephasing rate in a spinless Id system, we present 
a simple argument demonstrating the cancellation of the 
leading contribution 1/tcc to l/rj'. For the purpose of 



a qualitative discussion, we can consider inelastic scat- 
tering within the picture of a classical thermal noise 
with characteristic frequencies |a;| <C T. Let Uq{x) 
be a given realization of the fluctuating scalar poten- 
tial U^{x,t — 0) created by right-/lcft- movers at time 
t = 0. The peculiarity of the Id geometry is that in the 
clean case the initial density profiles of right-/left- movers, 
whose dispersion relation is linear, remain unchanged: 
they move with the plasmon velocity either to the left 
or to the right as a whole, and so does the fluctuating 
potential: 

U^ix,t) = U^{xTut) . (6.1) 

We now demonstrate that the phases acquired by the 
time-reversed waves propagating along a closed three- 
impurity loop are equal for arbitrary U^{x) and hence 
cancel each other in the interference contribution even 
before the averaging over fluctuations of U^{x). As we 
show in Appendix [Fl the velocity of electrons moving 
along the "Cooperon path" is renormalized by the in- 
teraction and coincides with u. Referring to the paths 
shown in Fig. [6^, we first follow the "forward" trajec- 
tory Xc — ^ Xa Xf, for a particle which starts 
moving to the right a,t t — and is backscattered at these 
points at times ta, tb, and ic, respectively, see Fig.[8l It is 
important to recall that the scattering amplitude 34 has 
been eliminated through the renormalization of the Fermi 
velocity. Therefore, on the right-moving parts of the for- 
ward trajectory (segments — > cca and Xb 0) the phase 
changes by (j)^ only due to the scattering off U~{x,t), 
while on the left-moving parts (segment Xa — > Xb) only 
the scattering off U~^{x,t) matters and yields the phase 
shift (j)^. Similarly, there are phase shifts for the back- 
ward path 0^ due to the scattering off U^{x, t). 

The phase acquired on the forward path between 
points Xa and Xb is given by 

(j)^= / U+(x)dx . (6.2) 

u J-uU 

Here we use tb ~ ta — tcjl. Electrons moving on the 
backward path interact with the same right-moving po- 
tential \J^(x, t) on the segments ^ Xb and Xa ^ 0, see 
Fig. [51 The phase (f>^ is written as 

1 

(f>^ = / U+{x)dx^(P^ . (6.3) 

^ J -2u{ta~tb+ta) 

Similarly, (f)^ — ((>+■ This is clearly seen in Figs. [5Ji and 
|9ji, where the regions of the potential profile that con- 
tribute to the phase shifts are identical at time tc for the 
forward and backward paths. As a result, the interaction- 
induced phase shifts cancel out in the total Cooperon 
phase {(j)^ + 4>-) — (0+ + <P-) for an arbitrary profile 
Uo{x). This means, in particular, that in a calculation 
of the conductance through a short piece of Luttinger 
liquid to third order in the impurity strength, which is 
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FIG. 8: Illustration of the quasiclassical electron dynamics in 
the presence of three impurities (whose positions are denoted 
by the dashed lines) and a right-moving component of the 
fluctuating scalar potential U'^{x,t). Top: the forward path 
corresponds to that in Fig. |6^. The position (black dots) and 
the direction of the velocity (arrows) of the particle, as well as 
the profile of the fluctuating potential, are shown at different 
times: (a) t ^ 0; (h) t = ta + 0; (c) t = h + 0; (d) t = tc + 0. 
The phase (f>^ of the electronic wave is accumulated when 
the particle moves to the left, that is between the snapshots 
(b) and (c). The segment of the fluctuating potential that 
contributes to the phase shift is marked in bold. 



the order at which the interference term arises, the lat- 
ter would remain unaffected by inelastic e-e collisions. In 
the calculation of the WL correction to the conductivity 
of an infinite quantum wire, we have, indeed, to include 
disorder in the interaction propagator, despite the bal- 
listic character of dephasing in the WL regime and in 
accordance with the argument made in Sec. IVI Al This 
is done systematically within the path-integral approach 
in the next subsection. 



D. WL dephasing in Id: Path-integral calculation 

The dephasing-induccd action S{tc, to) acquired by the 
Cooperon is accumulated on the classical (saddle-point) 
path, whose geometry for three impurities if fixed by 
two length scales, the total length of the path vptc and 
the distance between two rightmost impurities vpta < 
vptc/'^ (Fig. [To]). As shown in Appendix IB| the WL cor- 
rection can then be represented as 

AfJwl = -2(70 / dtc dtaP2{tc,ta)exp[~S{tc,ta)] , 

Jo Jo 

(6.4) 



FIG. 9: Same as in Fig.|8]but for the backward part, Fig.|6^. 
(a) t = 0; {h)t^tc-tt + 0; (c) t = tc - ta + 0; (d) t = U + 0. 
The phase (f)^ of the electronic wave is accumulated when the 
particle moves between the snapshots (a) and (b) and between 
the snapshots (c) and (d). 

where 

P2{tc, ta) = (1/8t2) exp(-ie/2r)e(ic - 2ta) (6.5) 

is the probability density of return to point x = after 
two reflections at points x = vpta and x — —VF(tc/2 — ta) 
(in this subsection, we ignore the difference between vp 
and u). The contribution Sij to the dephasing action 
associated with inelastic interaction between electrons 
propagating along the paths Xi{t) and Xj{t) is obtained 
similarly to higher dimensionalities 

X exp{iq[xi{ti) ~ Xj{t2)] ~ iuj^ti - t2)} , (6.6) 

where y^j,(a;, q) is the dynamically screened retarded in- 
teraction. Because of the HF cancellation of the bare in- 
teraction between electrons from the same chiral branch 
(f/4 processes), the screened interaction V{(jj,q) acquires 
the indices /i, v denoting the direction of motion of the in- 
teracting electrons, /x = sgni^, v = sgnxj. More specif- 
ically, the 04 processes have been incorporated in the 
renormalization of the bare Fermi velocity. It is worth 
stressing that if we would keep both g2 and (74 processes 
in V{u!, q), the dephasing action in Id could not be writ- 
ten in the form of Eq. (|6.6p . Indeed, in that case the 
dephasing action would generate the contribution of the 
diagrams of type (b) in Fig. 21 which is actually absent 
due to the HF cancellation (|5.8|) . 

The main steps in the derivation of Eq. (|6.6p are as fol- 
lows. First, the interaction is treated within the random- 
phase approximation (RPA) , neglecting closed electronic 
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FIG. 10: Illustration of electron dynamics governing the 
WL and dephasing: Time- reversed paths Xf{t) (solid) and 
Xbit) = Xf{tc ~ t) (dashed) on which the interaction-induced 
action S that yields dephasing of the Cooperon is accumu- 
lated. Dotted lines: the propagation of dynamically screened 
interaction. The interaction may change the direction of prop- 
agation upon scattering off disorder (as marked by a cross). 
Each interaction line gives a contribution to S proportional 
to {Nf — Nt)'^, where Nf^i, is the number of its intersections 
with the forward (/) and backward (6) paths. One sees that 
A^^ Nf, only due to impurity scattering in the interaction 
propagator. Interaction and electron lines lying on top of each 
other do not yield dephasing because of the HF cancellation. 



bubbles with more than two interaction legs, which is 
justified for the weak interaction, a ^ 1. Second, we 
perform an independent averaging of the RPA bubble 
and the electronic trajectories over the disorder. This is 
justified if the characteristic energy transfer uj ^ I/t^'' is 
much larger l/r. Further, we treat the thermal electro- 
magnetic fluctuations through which electrons interact 
with each other as a classical field, which is a valid ap- 
proximation provided that w is much smaller than T. As 
will be seen from the results of our calculation, the char- 
acteristic energy transfer uj ~ I/t^' does satisfy both 
these requirements for the considered case of weak inter- 
action, a <C 1. Finally, we will neglect the influence of 
the interaction on the velocity (u vs vp), which is again 
justified to the leading order in the interaction strength. 
Expanding V^^ to second order in a, we have 



■7) 



where the functions J-^u Fourier-transformed to (x, t) 
space are given by Eqs. (|C.2p . (|C.3p with 7 = l/r. The 
action Sij can then be written in a simple form: 

Sij^TTa'^vpT dti dt2Tf_i^[xi{ti)-Xj{t2),ti-t2] , 
Jo Jo 

(6.8) 

where, to first order in r ^, J-i^ii,{x, t) read 

T++{x,t) ~ 5{x + VFt){l~\t\/2T) , (6.9) 
T+-{x,t) ~ Q{vlt^ -x^)/4,VFT , (6.10) 

and ^__(x,i) = J-"++(— a;, i), T ^{x,t) = (a;,t). 

The total action is given by 5 = 2(5ff — S'fb), where / 



and h stand for "forward" and "backward" time-reversed 
paths (Fig.Iini). 

Let us first calculate S for the case when the disorder in 
interaction propagators (|6.9p . (|6.10p is neglected, r^^ = 
0. We get straightforwardly 

5ff = S'fb - ^c?Tt,l2 . (6.11) 

One sees that Ss reproduces the Golden Rule result, 
Eqs. (jgT^ . fOn)) . As we explain in Sec. lVIEi this rate is 
relevant for the AB dephasing in the multiple-connected 
(ring) geometry. However, in the present case of dephas- 
ing in a wire the self-energy processes (S'ff + Sbb) are 
exactly canceled in the action by the vertex corrections 
(S'fb + Sbf), yielding S' = 0. This confirms the conclusion 
of Sec. IVI CI for clean-system interaction propagators, 
the WL dephasing is absent. 

Therefore, the dephasing in Eq. (|6.4p is only due to 
the dressing of the dynamically screened interaction by 
impurities. A calculation of the dephasing action S to 
order Oir^^') is presented in Appendix [Cl with the result 



S{tc,ta)^2T:a^Tta {t,-2ta)/T . 



(6.12) 



The dephasing vanishes for ta — 0, tc/2, since in these 
cases the two paths / and b are identical. 

Substituting Eq. ((6T2| in Eq. JlH) we find for r^^ < r: 

1 f°° 1 
Actwi = -- (td / dtc exp{-tc/2T) 
4 Jo T 

/•*c/2 

X / dtaexp[-Wri,(t,/2-t,)/r] 
2 

r 1 



r^iirj In^--^ln(a^T), (6.13) 



where 



1/2 



In agreement with the above findings, l/r^' vanishes in 
the clean limit, -22. in contrast to the total e-e scattering 
rate, Eq. ((57TO)) . The logarithmic factor in Eq. (pl^ 
is related to the contribution of rare configurations of 
disorder in which two of the three impurities are anoma- 
lously close to each other {ta — > or tc/2): in these 
configurations the dephasing is strongly suppressed. As 
a result, the characteristic t in the integral (|6.13p are 
spread over the whole range between t'^^ and r. 

The scale Ti defined in Eq. (I6.14p marks the tem- 
perature below which the localization effects become 
strong. It is worth mentioning that the T dependence of 
ct(T) = <td(7') + Af7wi(T) starts to be dominated, with 
decreasing T, by the WL term rather than by (7d(T) at 
T ^ Ti/a, i.e., well above Ti for q <C 1. These results 
are illustrated in Fig. [11] 

Let us briefly mention what happens below Ti. In 
the strongly localized phase, the only peculiarity of Lut- 
tinger liquid is the T dependence of the scattering rate 
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a(T) 




strong 
localization 



1/x Ti=l/a^x Tj/a Ep 

FIG. 11: Schematic behavior of cr{T) on the log-log scale. 
Dotted line: the T-dependent Drude conductivity. Below 
Ti/a the WL correction, Eq. (|6.13|l . dominates dlna/dlnT. 
Below Ti the localization becomes strong. 

T^^(T) and the related T dependence of the localiza- 
tion length ^(T), persisting down to Tt{T) ~ 1 (see 
Sec. IIII BT|) : otherwise, the picture is similar to that in 
many-channel quantum wires or 2d systems. If a <^ 1, 
there exists a wide range of T where the localization 
is strong but the conductivity is a power-law function 
of T. We call this transport regime in a weakly disor- 
dered system, associated with interaction-induced diffu- 
sion over localized states, Power-Law Hopping (PLH); see 
Refs. 65,70.72,104,105 for details. An important point is 
that the lifetime of localized states due to inelastic e-e 
scattering can be calculated in the PLH regime at the 
Golden Rule level. At still lower T, the lifetime becomes 
infinite, so that the e-e interaction can no longer support 
hopping via localized states and the conductivity van- 
ishes at a certain critical Tc, see Refs. I65ll66l for details. 

E. Weak localization vs Aharonov-Bohm dephasing 

As shown in Sec. IVIDi the dephasing rate that cuts 
off the WL correction is much smaller than the Golden 
Rule result from Sec. |Vl which governs the damping of 
the single-particle Green's function, Eq. (|5.20p . One may 
ask if the Golden Rule rate I/tco is observable at all. 
The answer is yes: this type of dephasing makes sense 
in a clean Luttinger liquid in the ring geometry, where it 
governs the decay rate (t^^)~^ of Aharonov-Bohm (AB) 
oscillations4^i^°°i^°^ Indeed, the amplitude Gab of the 
first harmonic of the AB conductance is proportional to 
the product of the Green's functions in two arms of the 
interferometer: 

Gab {— dre'"'-^ 

I ^ Jo 

X 9+{Li,T)g+{-L2,-T)}^^^^g , (6.15) 

where Li 2 are the lengths of the arms of the ring. In the 
case of Li — L2 = L, the damping is solely due to the 
dephasing and is given by 

Gab oc exp(-L/uToo) = exp{-L/uT^^) , (6.16) 



i.e., = Tee. The difference from the WL correction 
is in the absence of the vertex parts of the dephasing ac- 
tion. Specifically, 5*12 = 5*21 = in the AB setup [where 
Sij is given by Eq. (j6.6p with i,j equal to 1 or 2 for tra- 
jectories in arms 1 or 2, correspondingly], because of the 
absence of interaction between the arms, see Fig. [121 On 
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FIG. 12: Illustration of electron dynamics governing the de- 
phasing of the first harmonics Gab of the AB oscillations in 
the ring geometry with two arms of equal length L. The posi- 
tion of a particle on the ring is parametrized by the angle (/)(t) : 
(l> — corresponds to the point at which the ring is connected 
to one lead, (p = ±7r — to the point at which it is connected to 
the other lead. The paths 0i(t) (solid line) and <^2(i) (dashed) 
traveling along arm 1 and arm 2, respectively, interfere at the 
second contact at time t = L/u. Dotted lines: the propaga- 
tion of the dynamically screened interaction. The interaction 
may change the direction of propagation upon scattering off 
the contact (as marked by a cross). Similarly to Fig. 1111 each 
interaction line gives a contribution to the AB-dephasing ac- 
tion proportional to (A'^i — N2Y , where is the number of 
its intersections with the trajectory in arm 1 or arm 2. One 
sees that for any dotted line the number of possible intersec- 
tions in the time domain Q <t < L/u is strictly one, yielding 
only the "self-energy" contribution to the action and no "ver- 
tex part" . Interaction and electron lines lying on top of each 
other do not contribute to the dephasing because of the HF 
cancellation. 



the other hand, the self-energy part Sn = S22 (analog of 
Sff and S'bb in the WL problem) is given by Eq. (|6.1ip 
with tc — L/u. An analo gous T^^ dependence of r^^ 
was also obtained in Ref. 11071 for a quantum wire cou- 
pled capacitively to metallic gates. The preexponential 
factor in Eq. (|6.15p depends on details of how the AB 
ring is connected to the leads. A particular AB setup — 
two Luttinger liquids weakly coupled to each other a t two 
points by tunneling — was thoroughly studied in Ref. IIO6I 
To summarize, the WL and AB dephasing rates in a Lut- 
tinger liquid are parametrically different in the spinless 
case studied in this paper; specifically, ^ t^'. In dif- 
fusive many-c hann el wires, an analogous difference was 
shown in Ref. llOSl 
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VII. 



FUNCTIONAL BOSONIZATION: 
ALL-IN-ONE APPROACH 



Equations (|6.13p . (|6.14p have been derived in two 
steps, emphasizing two key ingredients of the theory: 
the Luttinger-Hquid renormahzation of disorder and the 
Fermi-Hquid-Uke dephasing of electrons scattered by the 
renormahzed disorder. Although both effects are gen- 
erated by e-e interactions, they are of essentially dif- 
ferent origin: the renormahzation of disorder is due to 
virtual processes with energy transfers in the range be- 
tween T and A, whereas the dephasing stems from real 
processes with energy transfers smaller than T. Concep- 
tually, an advantage of the two-step approach is that it 
reveals the transparent physical picture of transport in 
the disordered Luttinger liquid. Alternatively, one can 
develop a unified approach where the virtual and real 
transitions are treated on an equal footing, which is what 
we do below. A fully bosonized approach, as well as a 
purely fermionic one, turns out to be not particularly 
advantageous for this purpose. Instead, in this paper 
we employ the method of so-called "functional bosoniza- 
tion" , adjusted to deal with the particulars of e-e scat- 
tering in one spatial dimension. The method was intro- 
duced for a clean Luttinger liquid by Fogedby^^ and by 
Lee and Cheui^ A transparent exposition and further 
development of the method in Id can be found in re- 
cent reviews by Yurkevich^^ and Lerner and Yurkevichi^^ 
For recent applications to single-impurity problems see 
Refs. ,109.110. . 

In essence, the functional bosonization is a technique 
based on the conventional Hubbard-Stratonovich trans- 
formation that decouples the density-density interac- 
tion term in the fermionic action. Bosonic fields ap- 
pear in this approach as auxiliary fields of the Hubbard- 
Stratonovich transformation, whose fluctuations medi- 
ate the dynamically screened interaction between elec- 
trons. The peculiarity of a clean Luttinger liquid is that 
the RPA approximation for the dynamically screened 
interaction is exactfiii This is a drastic simplification, 
so that the ease in which an arbitrary strong interac- 
tion in the clean case can be incorporated via the fully 
bosonized method is also the property of the functional- 
bosonization approach. In the disordered case, the 
functional-bosonization method has a technical advan- 
tage, especially useful for weak interactions, in that 
the conductivity is represented (in contrast to the fully 
bosonized approach) as a fermionic loop, albeit dressed 
by interactions in all possible ways. 



A. Functional bosonization: Background 

Our purpose here is to develop the functional- 
bosonization scheme in the presence of disorder and to 
apply it for the calculation of the weak-localization term 
(|6.13p in the conductivity. The quantity to begin with is 
the fermionic Green's function. In the absence of disor- 



der, we write it in the Matsubara representation for an 
electron moving to the right or to the left (-) as 



where 



iT 



1 



2vf sm]i[T:T {x / V p ± it)] 



(7.1) 



(7.2) 



is the free Green's function and the factor exp(— B^^) 
accounts for interactions between electrons. The func- 
tion Bfj_p, can be found by using the fact that, in the 
absence of backscattering due to e-e interaction, the 
Hubbard-Stratonovich transformation characterized by a 
decoupling field </?(a:, t) results in an identical representa- 
tion of g^{x^ t) as the average over (p{x, r) of a function 

9 " 

r 



5^(x,t)M^ 



3a^(2;,t) = (g^(a;,r)) 



(7.3) 



where 5^ (a;, r) describes the propagation of a particle of 
the same chirality in the field </j(x,t): 

[dr ~ it^vpdx - i(p{x,T)]g^{x,T) = 5{x)5(t) . (7.4) 



The angular brackets in Eq. (|7.3|) and throughout 
Sec. IVIII below denote the functional averaging over the 
field (p{x,t). The representation (|7.3p is equivalent to a 
local gauge transformation for the fermionic field of the 
right- /left-mover— 

ipf,{x,T) ^ iljf,{x,T)exp[ief,{x,T)] , (7.5) 

where the bosonic field 9^{x,t) is related to the decou- 
pling field ip{x,T) by 



{dr - i^J.VFdx)Of,{x,T) = ip{x,T) 



(7.6) 



The Gaussian character of the Luttinger-liquid theory in 
the clean case means that the averaging in Eq. (|7.3p is 
performed with an action in which all terms containing 
powers of >f{x, r) higher than second vanish to zero. The 
correlator {(p{x,T)(p{Q,0)) = V{x,t) then gives the dy- 
namically screened RPA interaction. 

As a result, the functions B^^{x,t) are represented as 
propagators of the fields 9f^{x, r): 

B^,{x,T) = {[e^{o,Q)-e^{x,T)]e,{Q,Q)) , (7.7) 

where for later use we also introduce the correlator of 
fields 6^ of different chirality, /Lt, = ±. It is convenient 
to put the bare coupling constant 174 = 0, shift the Fermi 
velocity accordingly (as discussed in Sec lIII A]) and deal 

with two different interaction propagators V++ ^ 

(see Appendix [X)) . Then the correlators i3^,y(a;,r) read 

B+±{x,r)^Tj2f^ (e^,.-.n„. _ 1^ 

n 

V+±{inn,q) 
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where f2„ = 27rnT is the bosonic frequency, B {x,t) = 

.Tjt), and B ^{x,t) = B^ {x,t). The integra- 
tion with Vfj,i, given by Eq. (jA.5|l [or, equivalently, by 
Eqs. (jA.7|) . (jA.8[) with a renormalized Fermi velocity, see 
discussion in Sec. IIII A] yields (see Appendix ID|) : 

B^,f,{x,T) = -ln77^(x,r) - y£(x,T) , (7.9) 



Bf,^^f,{x,T) = — ^C{x,t) 



(7.10) 



where 



ab = il-K f/2K, ar = {l- K^)/2K , (7.11) 
(tT/A)2 



and 

C{x, t) ~ In 



ih[7rT'(x/u + it) ] sinh[7rT(a;/u — it) 



(7.12) 



The functions 



sinh[7rr(a;/'yF ± ir)] 
u smh[7rJ {xlu± it)\ 

replace vp ^ u'va the free Green's functions g''^' {x, r) in 
Eq. dZI]). 



B. Functional bosonization: Diagrammatic 
technique 

The exactness of the RPA approximation in Id sug- 
gests the following diagrammatic technique. The Green's 
function (|7.ip for a right-mover propagating from point 
1 — (xi,Ti) to point 2 = {x2,T2) can be represented 
by the diagram in Fig. 1131 where the solid line is 

5+^(2^2 — xi,T2 ~ Ti) and the wavy line connecting the 
end points denotes the pairing of two fluctuating factors 
exp[i0+(2)] and exp[— i0+(l)] with the Gaussian fields 
61+ (2) and 61+ (1): 

(exp{i[0+(2)-^?+(l)]}> 

= exp [-B++(a;2 - xi,r2 - Ti)] . (7.14) 

The fields ^+(2) and ^+(1) enter the exponent of 
Eq. (|7.14p with different signs, according to the direction 
of the arrow on the solid line. 

Consider now a backscattering vertex which describes 
reflection of the right-mover from an impurity (Figs. [14] 
and 1151 before and after averaging over the fluctuations 
of respectively). The Green's function propagating 
from point 1 to point 3 upon backscattering on the static 
potential Ub{x) at point 2 is given by 

Sg{x3,Xl,T3 - Tl) ^ J dT2 J dX2 Ub{x2) 
X 5+^(2^2 - Xl,T2- Tl)g''°\x3 - X2,T3-T2) 



X ( exp{ ^[0+(2) - 9+(l) + 0_(3) - ^-(2) ] } ) 



-ie+(i) 



+19^(2) 



-B++(X2-Xi,X2-Xi) 



+ 




+ 



FIG. 13; (a) The Green's function of a right-mover prop- 
agating between points 1 — {xi,ti) and 2 = {x2,T2) in a 
given field ip{x,T) acquires two phase factors exp[i0+(2)] and 
exp[— iS+(l)] (denoted as wavy lines) at the end points. Solid 
line: the free Green's function g^\x2 — xi,r2 — ti). (b) The 
wavy line connecting the points 1 and 2 denotes averaging of 
the phase factors over fluctuations of 9+{l) and 9+{2). 



-ie+(i) 



i[e^(2)-e.(2)] 



ie_(3) 



+ 




FIG. 14: The Green's function with a backscattering vertex 
(marked by a cross) between the end points. The vertex is 
dressed by a local fluctuating factor exp{i[^+(2) — 6'_(2)]}. 



A direct consequence of Eq. (|7.5p is that the ver- 
tex at point 2 is dressed by the local fluctuating fac- 
tor exp{i[6'+(2) — 6'_(2)]}. The correlations of the fields 
at different points in x space in the last factor of 
Eq. (|7.15p account for e-e interaction. Apart from the 
averages (I7.14p of the fields of the same chirality, there 
appear averages of the type 

(exp{z[0+(2)-0_(l)]}) 

= exp[-B+_(a;2 - a::i,r2 - ti) - Mq ] , (7.16) 



where 



Afo = (02(o,O))-(0+(O,O)0_(O,O)) 



(7.17) 



(7.15) 



Note that the factors exp(±Afo) cancel out in any closed 
fcrmionic loop, so that in the diagrammatic technique for 
observables one can omit them from the very beginning. 
The Gaussian character of the correlations means that 
the exponent of the last factor in Eq. (|7.15[) is represented 
as a sum of the correlators B^^i, : 

( exp{ i[e+{2) - 0+(l) + 0_ (3) - 0- (2) ] } ) 
= exp [-S++(a;2 - a;i,r2 - n) 

+B+^{X2 - Xi,T2 - Ti) - B (X3 - X2,T3 - T2) 

+B+^{X3 - X2,T3 - T2) - B+^{X3 - Xl,T3 - Tl)] , 

(7.18) 

as illustrated in Fig. [13] [we used in Eq. (|7.18p the fact 
that B+_(2,2) = 0]. 
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B^3,l) 




6^2,1) 



B,.(3,2) 



B+.(2,2) 



FIG. 15; Illustration of the diagrammatic technique with the 
propagators B^ji,: the correction 5g{x3, xi,T3 — ti) [Eq. (|7.15|l ] 
to the Green's function due to the backscattering at point 2. 
The wavy lines denote averaging of the fluctuating factors 
exp(±6p) shown in Fig. 1141 The pairing of all the fields 9^ 
with each other according to Eq. (|7.7p yields the correspond- 
ing contributions _B^y to the factor (|7.18p . For ease of visu- 
alization, the ends of the free Green's functions are split off 
from the backscattering vertex at point 2. 



The appearance of the combination 0+{2) — 0-{2) for 
backscattering at point 2 in Eq. (|7.15p and the factoriza- 
tion of the correlations (|7.18p in the simplest case of a 
single backscattering vertex suggests a way to proceed to 
two or more vertices. In effect, the e-e interaction can be 
completely gauged out to the impurity vertices, as em- 
phasized in Ref. [osl . Specifically, one can, in principle, 
calculate observables (given by closed fermionic loops) 
perturbatively to any order in the strength of disorder 
and exactly in the strength of e-e interaction by attach- 
ing a factor 



exp{±i[0+{N)~0.{N)]} 



(7.19) 



to each backscattering vertex, where N labels the vertex, 
and pairing all the fields dfj.{N) with each other. The av- 
eraging over the fluctuations of 6^ {N) is then performed 
with the correlators (|7.9p - (|7.10p calculated for a homoge- 
neous system. The sign ± in Eq. (|7.19p should be chosen 
to correspond to the chirality with which an electron is 
incident on the impurity. 

The WL correction to the conductivity is given by the 
fermionic loop with 6 backscattering vertices (Fig. [5]), 
each of which should be dressed with the factors (|7.19p . 
However, taking disorder into account at this level would 
not be sufficient. As we showed in Sec. lVIl in order to cal- 
culate the effect of dephasing for spinless electrons, disor- 
der should also be included, to first order in the strength 
of disorder 7, in the loops that make up the RPA se- 
ries for the effective interaction. We recall that in the 
RPA approximation, substantiated in Sec. IVIl no more 



than two external interaction lines connect to each of the 
fermionic loops that constitute the interaction propaga- 
tor and the averaging of the interaction propagator over 
disorder is performed in such a way that no impurity 
lines connect different fermionic loops. As we show in 
Appendix [El in the language of the functional bosoniza- 
tion, the RPA approximation in the presence of disor- 
der translates into a disorder-induced renormalization of 
the correlators B^^{x,t). Specifically, the renormaliza- 
tion consists in adding a disorder-induced damping to the 
propagators V++ and V+_ in Eq. (|7.8p . The renormal- 
ized interaction propagators are derived in Appendix |^ 
[see Eqs. (|A.24l) ~ (jA.25p ] and are given in the Matsubara 
representation by 



V+-[iVLn,q) 



? 9 9 



= 92 



2TIVF q\L'^ +nl+\nn\{-t + 5-i) ' 

(7.20) 

q^vl + nl + \n^\[^ + 5^/2) 
+ +|f7„|(7 + 57) ' 

(7.21) 



where 7 denotes the transport scattering rate, Eq. (jA.lSp . 
and 57 is given by Eq. (|A.26p (at variance with the pre- 
ceding part of the paper, in this section we do not use 
the notation r for the transport scattering time in order 
to avoid a possible confusion with the Matsubara time) . 

Thus, our strategy in deriving Eq. (|6.13p via the func- 
tional bosonization is this: find the disorder-renormalizcd 
correlators B^,y{x,T) to order 0(7g|) and use them to 
average over fluctuations of the fields 0^ "attached" to 6 
impurity vertices of the 3-impurity Cooperon (Fig. [T6)) . 
The diagrams (b) and (c) in Fig. [5] are calculated in the 
same way. 



C. Functional bosonization with a 
disorder-renormalized RPA 

We now turn to the calculation of the effective 
(disorder-renormalized) propagators B^y{x,T) to order 
0(73!) by using Eq. (|7.8p and the interaction propaga- 
tors V++ and V+- given by Eqs. (fT^ . fT?^ . At second 
order in interaction, B^-\-(x,t) reads 



S4 



9l 



2 2 



27ruF 
- iqvpfl. 



2tt 



(-iO„ + qvpf {q^vj, + Ql + |f7„|7) 



1) 



(7.22) 



Let us represent as a sum = Bq + By + B^, 

where Bo{x) is the contribution of ri„ — 0, By{x, r) and 
S<(x, r) = S>(— X, — r) are the contributions of ri„ > 
and Qn < 0, respectively. In turn, By{x,T) is a sum of 
two terms: 



By=B 



dp 



Bt 



(7.23) 
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FIG. 16: The three-impurity Cooperon diagram with in- 
teraction effects encoded in the fluctuating factors exp(±0^) 
(denoted by the wavy lines) attached to the backscatter- 
ing vertices. Each impurity vertex A'^ = {xn,tm) is char- 
acterized by two fluctuating fields 9±{N), and similarly for 
iV = (xnjTn)- Pairing of all the fields with each other (simi- 
larly to Figs. I14I15P describes both the elastic renormalization 
of the three impurities and the dephasing of the Cooperon. 



where i?^' is the contribution of the double pole a,t qvp = 
if2„ and By is the contribution of the single poles at 

qVF = ±i{^l + ^nlY'^. 

The term coming from zero frequency does not depend 
on 7: 



2'KVl 



1 



2tt q^vj, 



,7rT|x| 

VF 



, (7.24) 



where a — g2/2T:vF, as before. The contribution of the 
double pole does not depend on 7 either, 



^^e(x) 

VF 



1 



(7.25) 



exp [2'kT(xIvf + jt) ] — 1 

and is an expansion of In 77+ (j7.13p in powers of a. Ex- 
panding the single-pole contribution to first order in 7, 
we have 

r, 27rT/A 
< m 

4 

7^1 



■In 



1 - exp[-27rT(|2:|/t;i=^ + ir)] 
1 



2vF l-exp[-27rT(|a;|/uF + iT)] 



7 



■F 



2nT 



VF 



where 



(7.26) 



(7.27) 



The sum (TTTf)) is defined for Rez > 0. At small \z\ < 1, 
it behaves as F{z) ~ zln|z|, which cancels the second 
term in Eq. (|7.26p at the ultraviolet cutoff r — > 0, 

Below, we will need the analytical continuation 
of B^i,{x,t) in the complex plane of r. In particular, 
an important contribution to the dephasing will come 
from the analytical continuation of F{z) to large negative 
Rez < 0. The latter can be done, e.g., by means of the 
integral representation 



F{z)= [ dz' ln[ 1 - cxp(-z') ] 
Jo 



(7.28) 



Piecing all the terms together, B++{x,t) with a 



disorder- induced correction of order 0(732) 
as 



is written 



B++ — B++\j=o + b , 
where i?+4_|^=o is given by Eq. (|7.9p . b 
7 I / 2'nTx 



bi{x,T 



In 
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IBttT 



\ VF 
1 



(7.29) 

= 61-1-62, and 

2ttT\x\ 

VF 



4sinh[7rT(a;/tii? + ir) ] sinh[7rT(.T/i'i? — zr) 



62(0;, t) 



7 



IGttT 



F 



2ttT 



VF 



27rT( ^ 

VF 



(7.30) 



Note that the second term in Eq. (|7.30p can be viewed 
as an x dependent renormahzation of a induced by dis- 
order, da = —a"f\x\/4:VF, as follows from a comparison 
with Eq. (|7.9p . This is consistent with the RG flow of the 
interaction constant However, Eqs. (fO(I)) . (|7^ 

contain more information and allow us to study the de- 
phasing missed in the RG equations. 

Consider now the disorder-renormalized propagator 
_B_l (x, r). To second order in g2, 



B+^{x,r)o,-Tj2 
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9l 



dq 

2^ nl + q^vj, 



AnvFq^vl + ni + \nr,hJ ' ^^'^^^ 

Expanding Eq. (|7.3ip to linear order in 7, we readily 

obtain a correction to B^ {x,t) [Eq. ()7.10p ] which is 

again given by Eqs. (|7.30p . (|7.30p but with the opposite 
sign, so that finally 

±\^^o±b. (7.32) 



B. 



B. 



D. Cooperon in a dynamic environment 

Now that we have found the propagators Bf^i,(x,T), 
let us turn to the calculation of the WL correction to the 
conductivity. The interaction induces the factor 



exp(-5c) = (exp[z(0/-0b)]) 



(7.33) 
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to the Cooperon loop, where 6f and 9i, are the phases 
accuniulated by an electron propagating along the "for- 
ward" and "backward" paths (Fig. [IT] represents the 
paths from Fig. [S^ in a way which is more convenient 
here) and the averaging is performed with the disorder- 
renormalized correlators ()7.32p . Again, Sc accounts for 
both the dephasing and the elastic renormalization of im- 
purities. The factor (|7.33p is the same for all 3 diagrams 
for AfTwi in Fig. [3 



+ B {XN ~ Xn',Tn - TN') 

- 2B+^{xN ~xn,,tm ~TN') , (7.39) 

and similarly for the pairings between N and N . The 
constant Mq given by Eq. (|7.17p cancels out in Sc ■ Note 

that on the closed contour {^^/^ = (^fo)' which is ac- 
counted for by the factor of 2 in Eq. (|7.36p . 



forward 




backward 



FIG. 17: Two time-reversed paths whose interference consti- 
tutes the three-impurity Cooperon. There are three backscat- 
terings off impurities on each of the paths (the positions of 
the impurities are marked by the dashed lines). The one- 
dimensional paths are stretched vertically for ease of visual- 
ization, with the "time axis" directed upwards for the forward 
path and downwards for the backward path. Note that this 
representation of the paths is identical to the graphs in the 
Berezinskii technique.^® 

The phases Of j, are given by 

Of - ^(-l)^+M^+(^)-^-(^)] , (7.34) 

N=l 



N=l 



where N,N are the impurity- vertex numbers, for the 
forward and backward paths, respectively; as shown in 
Fig. El The exponent in Eq. (I7.33P is represented as 



where 



Mff 



Sc = ~2{Mfj - Mjt) , 



(_l)A^+w'M(iV, N') - Mo , (7.37) 



M 



fb 



NN' 

l {Of Ob) 



. ^(-i)^+^Af (iV, N) - Mo , (7.38) 



NN 



and all the pairings between the phases associated with 
impurities N and N' are given by the combination 

M{N,N') = B++{XN-XN',TN-TN') 



E. Diffuson ladder and renormalization of disorder 

Before calculating the three-impurity Cooperon, it is 
instructive to look at the diffuson diagrams. More specif- 
ically, let us consider first the ladder diagrams without 
including self-energy corrections. 



1. One impurity line 

The simplest diagram for the density-density correlator 
of order 0(7) is shown in Fig. 1181 and written as 



di 



2 



l/T i-l/T 

dri / dfi 
^0 



g\'{xi,Ti)gl'{x -xi,T- n) exp[-M(l, 1)] 



X g^°\-x,,-f,)g'yix 



X, Ti - T) 



(7.40) 



The dashed line gives the bare strength of backscattering 
Wi?7o/2. In Eq. (|7.40p we only integrate over internal time 
variables, keeping other coordinates in (x,t) space fixed. 
The factor exp[— M(l, 1)], which is taken at the same 
spatial point, depends on ti — fi only: 



exphM(l,l)] = HA 



sm 

i 

' A 



^^ + A 



(7.41) 



(7.36) where we used the ultraviolet cutoff from Eq. ()D.7p and 



(7.42) 



Note that the interaction-induced phase factors cancel 
out at the external vertices which do not change chiral- 
ity. In particular, this means that in Eq. (|7.40p there is 
no renormalization of the velocity by the factors (|7.13p . 
The latter is a general property of fermionic loops with a 
larger number of impurities: the Green's functions con- 
necting to the external vertices with a small momentum 
transfer are characterized by the bare Fermi velocity vp. 
In the absence of e-e interaction, di is equal to 



4°) 



(7.43) 
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1/T 



FIG. 19: Analytical structure of the integrand in Eq. (|7.44|l . 
with a pole at ri = ix\/vF, for > and x < xi. Closing the 
contour of integration upwards, the integral along the real axis 
of the Matsubara time is represented as a sum of two integrals 
along the imaginary axis at ri = +0 and ri = 1 /T — 0. Using 
the periodicity with respect to the shift ri ^ ri ± 1/T, the 
sum gives a residue at the pole. 



FIG. 18: Density-density correlator of first order in the im- 
purity strength. The backscattering vertices are marked by 
the crosses. Pairing of the fluctuating factors exp(±0^) (de- 
noted by the wavy lines) with each other yields the factor 
exp[— M(l, 1)] in Eq. (|7.40|) . which describes the renormaliza- 
tion of the impurity strength by the interaction. 



where 



1/T 



/= / dTig''^\xi,Ti)g^^\x - xi,T - Ti) 
Jo 

_ T 1 

~ ~2v^ sinh{7rT[ {2xi - x)/vf + ir]) 

X [e(a;i)e(xi-a;)-e(-a;i)e(a;-a;i)] . (7.44) 

What is important to us here is that / is given by the 
pole on the classical trajectory ti = xi/vf of a right- 
moving electron (for xi > 0) or a left-moving hole (for 
xi < 0), where ti = —ivi (Fig. [TQ]). If we proceed to 
the O representation, the diagram for nonzero x and xi 
is written as 



1/T 



X [ e(xi)e(a;i - x)e{nn) + ei-xi)e{x - xi)ei-n„) ] 

and again is given by the (double) pole on the classical 
trajectory t — {2xi — x)/vp, where t — —ir and |2a;i — x\ 
is the total length of the trajectory from point x — 
to point X with backscattering at point xi. This simple 
consideration suggests that the two-particle correlators 
in the presence of a weak e-e interaction should still be 
determined by a close vicinity of the classical trajectory. 

To evaluate the integrals (|7.40p at ^ 0, we assume 
for definiteness that xi > and x < x\. The singularities 
of the integrand of Eq. (|7.40p in the upper half-plane of f\ 
for —0 < Refi < 1/T — are then a pole at fi = ixi/vp 
and a branch cut which comes from the factor (|7.4ip . The 
branch cut starts at fi = ti -\- i / h. and is sent upwards, 
as shown in Fig. [20k . The contour of integration over fi 
can be closed upwards, so that the integral over fi along 
the real axis from fi = to fi = 1/T is given exactly 



by the contributions of the pole and the branch cut. The 
contribution of the latter is proportional to ae and can 
be omitted in the limit of weak interaction, after which 
the remaining integral over ri reads 



11 (J_ 
2 \2vF 



(A/nT) 



2a„ 



dri 



sm 



sinh{7rr[ (2xi ~ x)/vp + ir]} 



sin{7rr[Ti — T + i{xi — x)/vpW 



(7.45) 



Here the contour of integration C in the complex plane 
of Ti runs anticlockwise around the branch cut along the 
imaginary axis of ti. Fig. [^Ob . For |cke| ^ 1 the sin- 
gularity in Eq. (|7.45p is almost a pole on the classical 
trajectory ti = xi/vp with ti = —iTi, so that, indeed, 
the system dynamics is determined by a close vicinity 
of the classical trajectory. We observe, however, that 
taking the interaction-dependent factor in the integrand 
straight at ti x\/vp yields an ultraviolet singularity 
that exactly cancels the factor (A/7rT)^"= in front of the 
integral in Eq. (|7.45p . which returns us to the noninter- 
acting result. It is thus the integration in a close vicinity 
t\ — x\/vf 1/2^ of the classical trajectory along the 
branch cut that gives the renormalization of disorder by 
interaction. If \ae\ ^ 1, Eq. (|7.45p reduces to 



di ~ (A/ttT)^ 



(7.46) 



The difference between Eq. (j7.46l) and the noninteracting 
result consists solely in the last factor which describes the 
renormalization of the disorder strength 70 — > 7, where 



7 = 70 (A/ttT) 



2Qe 



(7.47) 



2. "Two-impurity diffuson" 

Before generalizing to the Cooperon diagram, it is use- 
ful to consider a diagram of order 0(7^) with two par- 
allel impurity lines. Fig. [^T] What is new in this "two- 
impurity diffuson" as compared to Eq. (|7.40p is that now 
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ixJvt. 



Xi 
a) 



l/T 



b) 



FIG. 20: Analytical structure of the integrand in Eq. (|7.40|l . 
(a) The contour of integration over f\ is closed upwards where 
the singularities of the integrand for xi > and x < x\ 
are a pole at ri = ixi/vp and a branch cut which starts at 
ri = n + i/A and runs upwards, (b) The integral over t\ of 
the contribution of the pole in Fig. 120b is given exactly by 
the branch cut starting at ri = ix\/vF and running upwards 
along the imaginary axis of t\. 



there are factors exp(±M) which connect different impu- 
rities. Specifically, the interaction-induced factor is 

exp[-M(l,l) - M(2,2) 

-M(l, 2) - Af(l, 2) + Af(l, 2) -t- Af(l, 2) ] . (7.48) 

Each of the cross-terms (the last 4 terms above), when 
taken on the classical trajectory, may give a large contri- 
bution to the exponent of Eq. (I7.48p . An important ob- 
servation, however, is that on the classical trajectory they 
cancel each other exactly. This can be seen immediately 
since on the classical trajectory ti = fi and at this point 
M(l, 2) = M(i, 2) and M(i, 2) = M(l, 2). The cancella- 
tion can be verified for diffuson diagrams of higher order 
in 7 and is a manifestation of the mere fact that for the 
density-density correlator there is no dephasing in the 
diffuson channel. No dephasing is generated both for the 
bare correlators B^i,{x,t) and for those renormalized by 
disorder at the level of RPA. 



+ 




+ 



our purposes here to analyze the case of weak interaction. 
The analytical structure of the integral over four time 
variables, n, r2, n, and f2 can be simplified by selecting 
the singularities that are strongest in the limit \ae\ ^ 1- 
Let us assume for definiteness that xi > and X2 < xi. 
In the plane of ri we have for Imri > 0: a pole at ti = 
ixi/vp and three branch cuts starting at ti = fi + i/A, 
T2+i{xi—X2)/u, and f2+i{xi—X2)/u, as shown in Fig. [511 
For weak interaction, these branch cuts are "weak" in the 
sense that they are characterized by a small exponent ae ■ 
The approximation we make consists in neglecting all 
weak cuts and keeping only "strong" cuts characterized 
by exponents close to 1. Taking into account the weak 
cuts yields corrections small in ae ^ 1. 
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FIG. 22: Integration over the internal times in the two- 
impurity diffuson shown in Fig. [2T] (a) The analytical struc- 
ture of the integral in the complex plane of ti: the contour 
of integration runs around a pole and three vertical branch 
cuts, (b) The contour of integration in the complex plane of 
•fi runs around three vertical branch cuts. 

Picking up ri at the pole ri = ixi/vp yields the follow- 
ing singularities in the plane of fi for Imfi > 0: a strong 
branch cut at fi — ixi/vp and two weak branch cuts at 
fi = T2 + i{xi — X2)/u and f2 + i{xi — X2)/u. Neglect- 
ing the weak cuts and integrating along the strong one 
similarly to Eq. (|7.45p gives the renormalization of the 
impurity at point xi. During the integration along the 
strong cut, we pick up other factors in the integrand at 
the classical trajectory, i.e., at the starting point of the 
strong cut, which yields a cancellation of the cross-terms 
connecting two impurities. The remaining double inte- 
gral over T2 and f2 with the factor exp[— M(2, 2)] has the 
same structure as in the case of a single impurity and 
gives the renormalization of the impurity at point X2. 
The net result is that the two-impurity diffuson 3,2 for 
Ickel ^ 1 simply acquires two identical factors (A/7rr)^"= 
which renormalize the impurities: 



4") [K/nTf 



(7.49) 



FIG. 21: "Two-impurity" diffuson with the fluctuating fac- 
tors &qy[±6 f_i[N)\ (denoted by the wavy lines) attached to 
the backscattering vertices (marked by the crosses). 

As far as the integration around the classical trajectory 
of the two-impurity diffuson is concerned, it suffices for 



where 



4"' = 



{Tl2vFf 



sinh^{7rr [ {x + 2xi - 2x2)/vf + ir]} 
X [Q{xi)Q{xi - X2)Q(x - X2) 
+ Q{-xi)Q{x2~-xi)Q{x2~x)] . (7.50) 
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Self-energy corrections are incorporated in the same way 
leading to an additional factor d2 d2 exp(— 7|a;-|-2xi — 
2x2\/vf), where 7 is given by Eq. (|7.47p . This procedure 
is straightforwardly extended to the difFuson ladder with 
three or more impurity legs, so that the difference of the 
full diffuson ladder in the presence of e-e interaction from 
that in the absence of interaction is only in the change 
70-^7 [Eq. ^M]- 

F. Dephasing action and weak- localization 
correction 

Let us now return to the three-impurity Cooperon 
whose interaction-induced factor is given by 
Eqs. ((7351) . ((7:^ - ([7:^ . As the analysis of the 
two-impurity difFuson demonstrated, the analytical 
structure of the diagrams with increasing number of fac- 
tors exp(±A/) rapidly gets quite comphcated. The same 
example suggested, however, a systematic way to deal 
with the limit of weak interaction by ignoring in the first 
approximation all weak branch cuts. The integration 
around strong cuts then leads to the renormalization of 
impurities, while picking up nonsingular factors at the 
classical trajectory (the starting points of the strong 
cuts) generates the dephasing action. At this level, the 
calculation of the Cooperon is similar to that of the 
diffuson, with one essential difference: in contrast to 
the diffuson, the Cooperon factor exp(— S'c), Eq. (|7.33p . 
does not vanish on the classical trajectory and hence 
leads to dephasing. 

Another important issue is related to the definition of 
the classical trajectory for the Cooperon; more specifi- 
cally, to the question of what is the velocity that char- 
acterizes this trajectory. To answer this question, one 
has to carefully treat the factors r]± (I7.13P in the in- 
tegration over the internal times of the Cooperon. As 
shown in Appendix iFl the classical Cooperon trajectory 
is characterized by a single velocity and this is u. Having 
established the fact that there is the unique velocity on 
the whole classical trajectory of the Cooperon, we can, 
in the leading approximation, ignore below the difference 
between u and vp in the functions M{N,N'). 

Neglecting first the disorder-induced corrections 
b{x, t) in Eq. (|7.32p and summing up according to 
Eqs. (fr771) . (r05)) all 12 terms each of which is given by 
Eq. I|7.39p . we observe that the difference Mff — Mjb 
in this approximation "miraculously" cancels out to 
zero on the classical trajectory. This is in agreement 
with Eq. (|6.14p saying that the dephasing rate for the 
Cooperon vanishes in the limit 7 — > 0. 

Combining Eqs. (fr32| and (fr39| . the disorder- 
renormalized correlators M are represented as 

M^M\-y^Q + Ab, (7.51) 

where 6 = 61 -f 62 is given by Eqs. (|7.30p . (|7.30p . In order 
to evaluate the dephasing action on the classical trajec- 
tory, we analytically continue bi,2{x,T) in the complex 



plane of t. The analytical continuation of 62(2:, r) is ef- 
fected via Eq. (|7.28p . On the imaginary axis of r, the 
functions 61,2 at order 0(70^) read 

bi{x,it) ~ -^--f\x\ [{\x\ ~ VFt)Q{-\x\ + vpt) 

4 Vp 

+ {\x\+VFt)eHx\-VFt)] , (7.52) 
TT a^T 

b2ix,it) ~ ^-^-j [{\x\-VFtfO{-\x\+VFt) 

o Vp 

+ {\x\ + VFtf QhM - VFt)] (7.53) 

for both |x|/t)i? and \t\ ^ 1/T. Taking this limit is justi- 
fied for the calculation of the dephasing rate since upon 
averaging over the positions of impurities in Eq. (j6.4p 
the characteristic distance between impurities is vftJ^ or 
larger [in fact it is spread over the range between vft^^ 
and vft because of the logarithmic factor in Eq. (|6.13p ]. 
which in turn is much larger than vf/T according to 
Eq. dnUl). The sum of Eqs. ([73^ . ([735)1 gives 

, , TT a^T 2 2 2\ 

b(x,it) ~ -^^-f{x -vpt) 

o V p 

X [<d{-\x\+VFt) + <d{-\x\-VFt)\ . (7.54) 
Introducing a short-hand notation 

P{N,N') = ^—b[xN-XN',i{tN-tN')] , (7.55) 
7ra J 7 

where b is given by the asymptote (|7.54p . we see that 
/3(A^, N') = on the classical trajectory of the Cooperon 
for all TV, N' except for 

/3(1, 2) = /3(2, 3) = ~2Ut, ~ 2ta) . (7.56) 

Summing up the terms (I7.56P with the signs prescribed 
by Eqs. (|7.38p . (|7.5ip we obtain the dephasing action: 

5^"=^^ ~ 2Tia'T-fta{t, - 2ta) . (7.57) 

This result coincides with Eq. (|6.12p for the dephas- 
ing action derived via the path-integral calculation in 
Sec. IVIDI The integration along the cuts around the 
classical trajectory yields the renormalization of disorder 
on short spacial scales of order vf/T around impurities, 
similarly to the renormalization of the difFuson ladder in 
Eqs. (|7.46p . (|7.49p . The three- impurity Cooperon with 
both the renormalization of disorder and the dephasing 
included is thus given by 

C3 ~ (A/TTTf exp[-2na^T-ftaitc - 2ta)] , (7.58) 
where 

^(0) ^ ( loV {T/2VF? 

^ \2vfJ sinh^ [ttT {x^ / vf + iT)] 

X [e{xi)e{xi - X2)&ix3 - X2) 

X Q{x3 — x)Q{x3)Q{xi — x) 

+ Q{-xi)Q{x2 - xi)Q{x2 - X3) 

X e{x - X3)ei-X3)e{x - xi)] (7.59) 
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is the three-impurity Cooperon bubble without interac- 
tion, 

= 2(.Ti + X3 - X2) - X , (7.60) 

\xj:\ is the total length of the classical trajectory along 
which the density fluctuations propagate from point 
X = to point X. The length of the Cooperon loop 
vptc and the distance between two right-most impuri- 
ties vpta are related to the positions of impurities by 
vptc/'i = max{a:3,a;i} — X2 and vpta = \x3 — Xi\. Equa- 
tion (|7.58p describes the Cooperon bubble with scalar 
vertices — the current vertices are straightforwardly in- 
corporated below. 

To calculate the conductivity, one has to add also self- 
energy corrections to the Green's functions constituting 
the Cooperon bubble. In the presence of e-e interaction, 
this is done in essentially the same way as the renor- 
malization of impurities in Sec. IVII El and results in the 
additional factor exp{—j\x^\/vp), where 7 is the renor- 
malized rate of backscattering, which the Cooperon ac- 
quires on the length \x^\. The WL correction Actwi is 
then expressed in terms of C3 from Eq. (j7.58p as 

X y"dxiy"da;2y"dx3 cae-"^!^^!/"^! , (7.61) 

where the factor of 2 accounts for the contribution of 
the diagrams (b) and (c) in Fig. [SJ which give exactly 
the same contribution as the Cooperon, and the factor of 
(1/2)^ comes from the dressing of the current-vertices by 
the diffuson ladders, similarly to Eq. (jB.ip . The current 
vertices in Eq. (j7.6ip are simply given by ivp, since we 
formulate our model in terms of a linear electron spec- 
trum from the very beginning.-A-'^^ Integrating (|7.6H) we 
reproduce Eq. (|6.13p which was obtained in Sec. lVIDl in 
two distinct steps by treating the renormalization of dis- 
order and dephasing separately. Here, we have demon- 
strated how these two effects arise "side by side" from 
the technical point of view in the framework of the func- 
tional bosonization. Specifically, the analytical structure 
of the density-density correlator in the time domain is 
such that the dephasing comes from the classical tra- 
jectory, whereas the renormalization of impurities is as- 
sociated with the integration along strong branch cuts 
around the classical trajectory. 



VIII. SUMMARY 



We have developed two alternative approaches for a 
systematic analytic treatment of the problem. One is 
a two-step procedure which combines the bosonic-RG 
treatment of high-energy renormalization processes at 
the first step with the subsequent analysis of low-energy 
real processes in the fermionic language. The other ap- 
proach is based on the method of "functional bosoniza- 
tion" which makes it possible to treat both types of effects 
simultaneously. 

We have demonstrated that the notion of weak lo- 
calization is applicable to the strongly correlated one- 
dimensional electron system. This finding is of concep- 
tual importance, showing that the famous non-Fermi- 
liquid character of the Luttinger liquid does not pre- 
vent this system from exhibiting the features character- 
istic of conventional mesoscopic electron structures. Our 
approach thus provides a framework for systematically 
studying the mesoscopic phenomena in strongly corre- 
lated electron systems. Further directions of research 
include mesoscopic fluctuations and noise in a disor- 
dered Luttinger liquid, spin-related effects (in particular, 
the magnetoresistance), several-channel quantum wires, 
transport in a disordered system in the limit of strong 
interaction, transport in a "granular Luttinger liquid" 
(strong impurities). 
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In this paper, we have studied the transport proper- 
ties of interacting spinless electrons in a disordered quan- 
tum wire within the framework of the Luttinger-liquid 
model. Our main result is the weak-localization correc- 
tion (|6.13p . governed by the dephasing rate (|6.14p . 



APPENDIX A: POLARIZATION OPERATOR 
AND SCREENING IN DISORDERED WIRES 

In this Appendix, we derive the polarization opera- 
tor and the dynamically screened RPA-interaction for a 
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disordered Id electron system. As everywhere in the pa- 
per, we consider spinless electrons. The Luttinger-liquid 
power-law singularities only show up in the renormalized 
strength of disorder 7 and do not make any other differ- 
ence in the derivation below. 

We begin by considering the clean case. According 
to Dzyaloshinskii-Larkin theoremjiii the RPA equations 
for screened interaction are then exact. The propagators 
of interaction between right-movers (1^++) and between 
right- and left-movers {V-^ ) obey 



y++ = .94 - .94 n+ V++ - 52 n_ 
v+- = .92 - .92 n+ V++ - 54 n_ v+_ 



(A.l) 
(A.2) 



where 



1 



271 qvp ^ iflr, 



are the chiral polarization operators for right(-|-)/left(-) 
movers, fin = 2'KnT is the bosonic Matsubara frequency. 
The total polarization operator is a sum of the chiral 
terms: 



n(irj„, q) 



2 2 

q vp 
P^r~2 — 



(A.3) 



where p = l/irvp is the thermodynamic density of states, 
so that n(0, q) = p. For later use we also introduce the 
two-particle propagator 



D±{iVtn,q) 



irp 



which satisfies 



dp 
2^ 



\Qn\ i^n T qVp 



(A.4) 



G±''(iem + + q)G^± {i£m,p) = D±(irj„, q) 



X [e(-e„)e(e™ + f7„) + e(e„)e(-e„, -f7„)] , 

where G'£\iem,p) = (i^m T vpp)~^ are the free Green's 
functions for a given chirality and the ultraviolet momen- 
tum cutoff is sent to infinity. The polarization operator 
and the two-particle propagator are related to each other 
via 



n±(zf7„,g) = | 



27r 



The solution to Eqs. (jA.ip and (jA.2[) for arbitrary (72 
and 54 reads?^ 



54 



52^) 



1 



qvp 



(.94 ilZJ n •(-! I 

2'KVF i\ln + qvp 



V, = 



X {nl + q\l)i{nl + q^u^) , 

= g2{nl + q'vl)/{nl + q^u') 



where 



Vp 



(1 -I- gi/^irvp) - {g2/2TTVpY 



1/2 



For (74 = g2 we have 



y++ = y+- = 92 



ni + u^q^ 



u — Vp. 1 



92 
TTVp 



(A.5) 
(A.6) 



while for g^ ^ the propagators V++ and are dif- 
ferent from each other: 



gl qvp{iVtn- qvp) 



++ 



y+- = 92 



2Trvp Ql + vP'q^ 

nl + q^vl 



ni + u^q^ 



u = Vp 



1 ~ {g2/2TTVpy 



1/2 



(A.7) 
(A.8) 
(A.9) 



The interaction between left-movers is given by 1/++ with 
9 -9- 

Let us now turn to the disordered case. We only 
take into account the backward scattering induced by 
disorder, since the forward scattering can be completely 
gauged out in the calculation of the conductivity. We av- 
erage the RPA bubbles over disorder separately by using 
the ladder approximation, see Sec. IVI Dl The backscat- 
tering off disorder makes it necessary to introduce two 
indices p,v = ± in the polarization bubble Xl/^jy, corre- 
sponding to the chirality of the vertices of the bubble. 
The RPA equations for the interaction propagators in 
the presence of disorder read 



y++ = 94 - (.94n++ + g2U^+)V+ 
- (52n__-f g4n+_)F+_ , 



(A.IO) 



v+- = .92 - (.92n++ + 5i4n_+)y+_ 

- (54n__ + <?2n+_)y++ , (A.ii) 



where 



Il^^{inn,q) 



2tt 



■Df_,^{inn,q). 



We now calculate the two-particle propagators Dfj^i, in 
the RPA. Performing the analytical continuation to real 
frequencies ifln uj + iO, we introduce the retarded two- 
particle propagator with only self-energy 
impurity lines included. 



{DU-,q)r 



mp 



to — qvp + i"f/2 



where 



7 = 1/t 



(A.12) 



(A.13) 



and r is the transport scattering time. Note the factor of 
1 /2 in front of 7 in Eq. (|A.12[) . which reflects the fact that 
only backscattering off impurities is considered, i.e., — 
2t, where is the quantum scattering time, in contrast 
to the case of isotropic scattering, where Tq — t. The 
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propagator [D^_{u;, (z)]*-"-' is given by Eq. (|A.12[) with the 
change q —^ ~q, whereas the nondiagonal propagators 

The equations for dressed two-particle retarded prop- 
agators read 



D^_{uj,q) = [D^_{oJ,q)] 



(0) 



D^_{uj,q) = [D^^{^,q)f^^^D^^ico,q) 



and 

D^+iu;,q) = D^_{oj,-q) , 

Solving Eqs. (|A.14p - (|A.17p . we get 

UJ — q + 17/2 



Dj_{LU,q) = inp- 



Df_{LO,q) = -Tip 



7 L 

2 Cj2 _ q2y2^ 



lUJ"f 



(A.14) 
(A.15) 



(A.16) 
(A.17) 



(A.18) 
(A.19) 



and 



n__(c^,g) 

n+-(c^,g) 

n(c^,g) 



n++(c^,-g) 

P q'^Vp — qvpuj — 
2 



q'^vj. 



U^+{uJ,q) 
p 



2 q Vp — LO^ — iujj 



? 9 
q Vp 



^-^ q^'vi, — UJ 



(A.20) 
(A.21) 

(A.22) 



Equation (|A.22p reduces to Eq. (|A.3p in the clean limit 
7—^0 and acquires the conventional diffusive form 

W^'^icj, q) = pDq^/{Dq^ - iuj) 

in the limit lu <ti ^,qvp , where the fd diffusion constant 
D = vlh. 

From Eqs. (|XTO)) . ([XTT]) . (|X20|, and ((X2T|) . we ob- 
tain the retarded interaction propagators for 52 = <?4: 



92 



i + 52n 



9 2 2 
9 — W 

32-57^ 



ZW7 



with u given by Eq. (|A.6|) . For 94 = we have 



(A.23) 



y^^ - -9l 



n 



9l 



(f + .g2n+_)2-.g|n__n+_ 

q^v^p — qvpuj — 



2tivp q^u^ — ixP- — icij(7 + ^7) 



(A.24) 



V^^ = 92 



92 



1 + .92n+- 



(i + g2n+_)2-52n__n++ 

<^v\ — — iw(7 -I- ^7/2) 



cp'u^ — a;2 — iw(7 -|- 6^) 
where u is given by Eq. (|A.9p and 

92 



-7 



27rz;i? 



(A.25) 



(A.26) 



APPENDIX B: WEAK-LOCALIZATION 
CORRECTION TO THE CONDUCTIVITY IN ID 

In this Appendix, we calculate the WL correction to 
the conductivity in a noninteracting Id system, assuming 
that the phase coherence is broken by an external source 
of dephasing and the dephasing time is small; specifically, 
T0 <^ T. Under this condition, the conductivity can be 
expanded in powers of t^/t. The lowest-order diagrams 
are depicted in Fig. The sum of the diagrams with two 
impurity lines, shown in Fig. [71 yields a zero contribution 
to order (ei^r)"^. 

We formally consider only backward scattering, i.e., 
the impurity line corresponds to the correlation function 
of the backscattering random potential 

L-i(|%(g)|2) = l/2^pr 

(L is the system size) and changes chirality at every im- 
purity vertices. For the anisotropic scattering, we have 
T — Tq/2, where r is the transport time and Tq is the 
quantum lifetime. The retarded single-particle Green's 
function expressed in terms of t reads 



(e,p) = {eTvpp + i/2Tq 



{e^ypp + i/Ar) ^ 



The anisotropic scattering also introduces the current 
vertex renormalization: each current vertex is multiplied 
by a factor r/r^ = 1/2. 

The sum of the diagrams (b) and (c) in Fig. [5] with an 
impurity line covering the two-impurity Cooperon (two 
crossed impurity lines) is equal to the contribution acs 
of the diagram with the three-impurity Cooperon. The 
total three-impurity WL correction is thus given by 

ACTwI = 2 (7(73 • 



The expression for <Jca reads 



= -2 ( - 



2 2 

e Vp 



(B.l) 

where is the Fermi distribution function, the over- 
all minus sign reflects the fact that the product of the 
current vertices is negative in the Cooperon with three 
backscattering lines, and the factor of 2 is related to two 
possibilities of setting chiralities at the vertices: (-1 — ) 
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and ( — h). The function J{Q) describes the fermionic 
loop: 



J{Q) = 



1 [ ^ [ ^ f 

(27rpT)3 J 2tt J 2tt J 2tt 

X G^ip)G'^ip,)G^ip,)G^i-p + Q) 

X Gi{p)G^{-p2 + Q)Gii-pi + Q)G^{-p + Q) 



(B.2) 

Using the equahty 

Gl{p)Gi{p)^2^r[Gl{p)-Gi{J>)], 
we simphfy Eq. (jB.2p : 

J{Q)=AT^\P{Q)\'[P{Q)+P*m , 

where 

Now we introduce a phenomenological dephasing time 
through an additional decay of the Green's function: 

(e,p) ^ (e T vfQ + ij^T + il2T^Y^ , 

which yields 

1 1 



PiQ) 



2t wfQ + 1/2t + 1/r^ 



(B.3) 



1 



'^^^^ t\2t^ [{vfQY + [1/2t + 1/t^YY ■ 

(B.4) 

Substituting Eq. (|B.4|) in Eq. (|B.1|) . we find for T0 < t: 

1 



Ao- 



wl 



(B.5) 



where ctd = e^vprj'K is the Drude conductivity. 

Equation (jB.Sp for the weak-localization correction can 
be cast in the form of an integral over the time needed 
for a particle to return to the starting point after two 
backscatterings: 



r-OO 

AfTwi = -2o-D / (iic^2(ic)exp(-tc/T0) 
Jo 



(B.6) 



Here 



i < T 



is the return probability after two reflections (hence the 
subscript 2) for a particle which starts moving cither to 
the right or to the left and the overall factor of 2 in 
Eq. (|B.6p comes from the summation over the sign of 
the velocity at the starting point. 

Let us now turn from the phenomenological factor 
exp(— tc/T^) to the actual dephasing factor. The latter 



depends not only on the total time but also on the ge- 
ometry of the Cooperon path. We thus have to perform 
averaging over the geometry of the path at given tc- For 
this purpose we need the probability density of return 
P2(tc,ta) for a particle which starts moving to the right 
at point X — under the condition that the total length 
is vptc and the first reflection occurs at point x = vpta- 
This function is written as 

/•OO 

P2{tc,ta)^VF / dhPl{ta)Pl{tb-ta)Q{tb-ta) 

Jo 

X exp[-{t^-h)/2T]Q{tc-tb)S[x{Q] . (B.7) 



Here 



is the probability density of being reflected at point vpt 
(again, note that r is the transport scattering time and 
only the backscattering is taken into account). The factor 
exp [—{tc — tb)/2r ] in Eq. (|B.7p describes the probabihty 
of avoiding the backscattering on the last segment of the 
path and 

x(tc) = 2vF{ta - 4) + Vptc 

for tc > tb- We obtain 

P2(tc,ta) = ^e-W2- 



X / dtbeitb-ta)e{tc~tb)d[2{ta-tb)+tc] 
JQ 

= A Oitc ~ 2ta) . 

For the dephasing action which depends on t and ta we 
have 

/"OO /^OO 

AfTwl = -So-D / dtc / dtaP2{tc,ta)exp[-S{tc,ta)] , 
JO JO 

which reduces to Eq. (|B.6[) for the phenomenological ac- 
tion S{tc,ta) = tc/r^. 



APPENDIX C: DEPHASING ACTION IN THE 
PATH-INTEGRAL APPROACH 

In this Appendix, we present details of the calculation 
of the dephasing action defined by Eq. (|6.6p and given in 
the final form by Eq. (|6.12p . For the imaginary parts of 
the retarded interaction propagators (jA.24p . (|A.25p which 
enter Eq. (|6.6p . we have, retaining only the terms of order 

Attvf [q Vp ~ ^ j + 7 



Next, we introduce functions Tj^±{tj}^q) which obey 



ImFf±(w,g) 



9l 



ujT+±{u,q) 



(C.l) 
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Putting 32 = in and doing the Fourier transform 
to (x, t) space, we get 



^ p-7|t|/2 



7 



16u 



(C.2) 



(C.3) 



Here x± = x ± f pi, 

/++(^) 



^/i(2Vi)e(^) , 

/o(2z)e(z) , 



/o(2:) and Ii{z) are the modified Bessel functions. Since 
the dephasing action becomes of order unity at times 

t ^ T^'' <C r, we expand J>+(x,t) and to 

the first order in the disorder strength 7, which yields 
Eqs. (lOll. dCT) . 

Let us now turn to the calculation of the dephasing 
action S{t,ta) which can be represented as 

S{t, to) = S'ff + S'bb — -Sfb + -Sbf = 2(5ff — S'fb) , 



where Sij are given by Eq. (|6.6p . We split the forward 
path into three segments corresponding to the right- (R) 
and left- (L) moving parts, where the coordinate Xf(t) 
behaves as 

If, R: Q<t<ta] Xf{t) = VFt 

Ilf, L: ta<t<ta+*-^; Xf{t)=VF{2ta-t) 

Illf, R: ta + ^<t<tc; Xf{t) = VF{t-tc) . 

Similarly, for the backward (time-reversed) path we have 
Xb{t) = Xf{tc - t): 



lb, 
lib, 

Illb 



< i < 



Xb{t) 



-Vpt 



R 



Xb{t) = VFi2ta-tc + t) 
ta <t < tc; Xbit) = Vpitc - t) . 



The "self-energy" part Ss of the action is written as 

3 6 



5ff 



TJ« , Jff = ^ J„ + 2 ^ J„ , (C.4) 



71=1 



n=4 



where each of the terms Jn corresponds to the integra- 
tion over a certain pair of the forward-path segments (for 
brevity, we set w^? = 1 in the arguments of below): 

If — If : Ji — vp dti / dt2 
Jo Jo 



X T++{ti - t2,ti - t2) , 

Ilf - Ilf : J2 ^ VF J dti J dt2 

X T^-{-h+t2,h~t2) , 

Illf - Illf : ^ VF dti / dt2 

X J^++{ti - t2,ti ~ t2) , 

If - Ilf : 2 J4 ^ VF dti / dt2 

X [J^+^{ti+t2-2ta,ti-t2) 

+ T-+{-ti-t2+2ta,ti-t2)] , 

If - Illf : 2 J5 ^ vf dti / dt2 

X [T++iti-t2+t,,ti^t2) 

+ -^++(^^1 + ^2 — ^c, " ^2) ] , 

Ilf - Illf : 2 Jg ^ vf dti / dt2 

Jta Jta+ta/2 

X [T-+{-ti~t2+2ta+t,,ti-t2) 

+ T+-{ti+t2-2ta~t,M~t2)] , 

The "vertex" part 5fb of the action is given by 



S'fb — -T^^TJfh , Jfb — Jn 



(C.5) 



n=7 



where 



/■ta rtj2~ta 

If - lb : J7 vf dti / dt2 
Jo Jo 

X J'+-{ti+t2,ti-t2) , 

If - lib : Js ^ vf dti / dt2 

Jo Jta/2-ta 

X j:-^^{ti-t2~2ta+t^,ti-t2) 
nta rtc 

If - lllb : Jg = vf dti / dt2 

Jo Jto-ta 

X T+-{ti +t2 -tc,ti ~t2) , 

pta+tc/2 pta/2-ta 

Ilf - lb : Jio = VF dti J dt2 

X T^^{-ti+t2+2ta,ti-t2) , 

pta+ta/2 nta-ta 

Ilf — lib : ^11 = vp dti / dt2 

Jta Jta/2-ta 

X J^-+{~ti~t2+tc,ti~t2) , 

/ta+tc/2 pta/2-ta 
dti / dt2 

X T^^{-ti+t2+2ta,ti-t2) , 

pic — ta 

Illf - lb : Ji3 ^ vp dti / dt2 

Jta/2 + ta Jo 

X J^+-{ti +t2 -tc,ti -t2) , 
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Illf — lib '■ Jii ~ vp I dti / dt2 

Jta/2+ta Jta/2-ta 

X :F++{h~t2-2ta,h-t2) 
Hit - Illb : Jl5 = VF I dti / dt2 

Jta+ta/2 Jta-ta 

X T+-{ti +<2 -2tc.,ti -t2) 



Carrying out the integration yields: 
Ji 



and 



J2 
J3 

J9 



ta 

2 ' 

k 

4 ' 

1 /tr 



2 V 2 

i 

7- 



^8 I2 "^'^ 



1 Ac , \ Ac , 

2 2 "'"j^^T 2 



Jll - ^ 

J12 

Jl4 
J5 



tc 

— — f 

2 ^ 



"2" ^"^7 V 2" ~ ' 



4 I 2 " 



1 Ac , \ ta ft, 

2 2 "'-^ 2 



= J7 = Jl5 = , 



t 

Jff = -+1 — 

^ 2 ^ ^ 8 ' 



Ji 



fb 



t_c 2 

2 4 



8ta 



Js — JO) ~ '^ita ( — ta 



tc 



(C.6) 

(C.7) 
(C.8) 



We see that the difference Js — Jfh vanishes at 7 ^ 0, i.e., 
in the absence of disorder in the interaction propagators. 
Finally, we arrive at the expression for the total action 



S[t,ta) = 2X 



9l 



T{Js - Jfb) 



al rr . (tc . 



(C.9) 



Using a = g2/2TrvF, we obtain Eq. (|6.12p . 



APPENDIX D: FUNCTIONAL BOSONIZATION 
IN THE CLEAN CASE 

In this Appendix, we calculate the correlators 
Bfi^{x,T) in a clean Luttinger liquid for an arbitrary 
strength of interaction. Following the discussion in 
Sec. nil Al we consider only the §2 interaction, while the 
g4 interaction is accounted for in the shift of the Fermi 
velocity, Eq. (|3.12p . Substituting the interaction prop- 
agator V++{ifln,q) given by Eq. (|A.7[) [with vp under- 
stood here and below as the renormalized velocity (|3.12p ] 
in Eq. ([7^ . we have 



B..(x,.,.-ri:/| 



qvp 



2-KVp ni + q^v? 



piqx — iQ,nT _ 1 

X = 60 + S> + 6< , 



-iVln + qvp 



(D.l) 



where ;Bo(x), By{x,T), and B^{x,t) are the contributions 
of terms with f2„ — 0, ri„ > 0, and 0„ < in the 
Matsubara sum over the bosonic frequency r2„ — 27rnT, 
respectively. The sums over r2„ > and il„ < are 
related to each other as 



S<(a;,r)=S>(- 



-r) 



(D.2) 



To accurately treat the ultraviolet cutoff for an arbi- 
trary ratio u/vp, we introduce a finite range of interac- 
tion d, so that 52(9) depends on the transferred momen- 
tum q and vanishes for qd ^ \. The RPA approxima- 
tion, leading to Eq. (jP.ip , is exact in the Luttinger-liquid 
model for arbitrary (72(9). In the integrand of Eq. (|D.1|1 . 
52(9), vpiq), and u(q) are related to each other at given 
q in precisely the same way, Eq. (|4.2p . as in the case 
of 52 (9) — const. Everywhere below in this Appendix, 
except for Eq. (|D.6p . when writing 52, vp, or u, we un- 
derstand them as taken at g = 0. 

The term in Eq. (jP.ip with il„ = at > d reads 



Bo{x) 



2-kT . 



2u2 



Vp 



(D.3) 



In the integration over q at > 0, there are poles at 
q = iiln/vp and q = ±ifln/u. The first only contributes 
to the integral if f2„ < vp /d, whereas the latter — only if 
fin ^ u/d, which sets the limits of summation over n: 



g-J2„/Ajg-0„(|a;|/M+ir) _ J 



71=1 



(D.4) 



where there appear two ultraviolet scales Ay — vp/d and 
A = u/d and we approximate the cutoffs by the expo- 
nential factors exp(— r2„/Au) and exp(— fJ„/A) (the exact 
shape of the cutoffs plays no role in the infrared physics 
we deal with in this paper). 
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Using for ^ 1 and Re z ^ 1 /N the formula 

N 



1 



-ln(l-e-^)-lnA^ , 



we write By{x, r) for ^ d as 

A 1 - exp [-27rT(a;/u + ir)] 



B>{x,t) = e(a;)ln 

vp — u 
2u 



A„ 1 - exp [-2ttT{x/vf + ir)] 



In 



27rT/A 



1 -exp[-27rT(|a;|/u + iT)] 



(D.5) 



Combining Eqs. (|D.3|) . (|D.5|) . and (|D.2p . we obtain 
-B++(a;,r) given by Eq. (|7.9p . Note that the factor vp/u 
in 77+(a;, r) [Eq. (|7.13p ] comes from the ratio A/A^ in 
Eq. (jP.Sp . The remaining cutoff in C{x, t) is given by A. 

The Green's function [Eq. (|7.ip ] of a right-mover then 
reads for \x\ ^ d: 



9+{x, t) = .gf ^ (x, t) exp [-B++(a;, r) ] 
iT 1 



2m sinh [■nT{x/u + ir) ] 

(ttT/A)"" 



{sinh [ttT^x/u + it) ] sinh [7rT(x/u — ir) 

Note that the Fermi velocity vf drops out completely in 
the product g'^^\x,T)rij^{x^T). The Green's function in 
the real-time domain, Eq. (|5.16p . is obtained by means 
of the Wick rotation t ^ it + sgnt/A. 

The calculation of {x^t) is performed in a similar 

way: substituting Eq. (jA.Sp in Eq. (|7.8p . we get 



B+-{x,t) = -T 



dq e^q^-ii^^r _ 



1 



(D.6) 



Note that the denominator fl^ + q'^Vp in Eq. (|7.8p is 
canceled by the numerator in Eq. (jA.8|) . so that vp dis- 
appears from the problem right away and the integration 
over q in Eq. (|D.6|) is determined at ^ solely by the 
plasmon poles q = iifln/u. For |a;| ^ d, we arrive at 
Eq. (jTTOll . 

In Sec. lVIII and in Appendix [El we do contour integrals 
in the complex plane of the Matsubara time r and should 
be careful about the behavior of Bf^^(x,T) at a; ^ for 
T 0. When integrating in the complex plane of r, it 
is convenient to have the ultraviolet cutoffs in Bf^,j{x,T) 
that are independent of t. Our choice is to use Eq. (|D.4p . 
which yields a cutoff which depends on x: 



A 



sgn x = X + d sgn x 



(D.7) 



i.e., |a;| — > |a;| -I- d in the propagators B^^{x,t) in 
Eqs. (|7.9p ~ (|7.10p . and similarly in the disorder- induced 
corrections to them in Eqs. (|7.30p . (|7.30p . 



APPENDIX E: "DIRTY RPA" IN THE 
FUNCTIONAL BOSONIZATION 

In this Appendix, we demonstrate that the cor- 
relators (|7.7p in the presence of disorder are mod- 
ified according to Eq. (|7.8p with the effective RPA 
interaction V^,y (|7.20p . (|7.2ip whose dynamical proper- 
ties reflect the disorder-induced backscattering ("dirty 
RPA"). For this purpose, we directly calculate the 
disorder-induced correction to the averages of the type 
(exp{i[0^(a;, r) — 6,^(0, 0)]}) in the framework of the di- 
agrammatic technique formulated in Sec. IVII Bl In the 
expansion in powers of 7, the rules of the technique pre- 
scribe that impurity vertices are connected with each 
other via the correlators Bfj_^{x,T)\^^o in a spatially ho- 
mogeneous system. To obey Wick's theorem, multiple 
fermionic loops can only be present in this technique if 
each of them contains backscattering vertices, since all 
fermionic loops without impurities have been accounted 
for by the "clean RPA" propagators. We show below that 
the result for i?^^ obtained this way coincides for a <C 1 
with the result of the "dirty RPA", Eq. (|73^ . 




(0,0) 



(X,T) 



FIG. 23: Diagram in the functional-bosonization technique, 
which describes the addition of disorder in the RPA interac- 
tion propagator. Two backscattering vertices (0, 0) and {x, t) 
are connected to each other via the fermionic bubble with a 
single impurity at point x\. 

Let us consider the average of first order in 7, 



,i[e+(x,r)-e+(o,o)] 



dxil{x, T,Xl) , 



Ae+{x,r}-e+{o,o)] 



7=0 



(E.l) 



given by the diagram in Fig. [531 We are interested in 
the pairings between the fields 9+ (0,0) and 9+{x, r) that 
connect up the points (0,0) and (x,r) via the fermionic 
bubble that contains two backscatterings at point xi: 



/(x,r,xr) = ^e---(-) 



l/T 



dri 



l/T 



dT[ 



35 



x.gf (0,r(-n)/°)(0,n-rO 

^ g-B++(0,Ti-T{)-B__(0,Ti-r0+2S+_(0.Ti-r0 

,-B++{xi,t[)+B^^(xi,t[)+B-)- + {xi,ti)-B-)-^{xi,ti) _ 



^B+^{x-xi,r-T[)-B-f-^(x-xi,T~T[) 



X e 



-i3++(x— xi,T— Ti)+_B+_(a:— xi,T— Ti) 



(E.2) 



i\x,\/u 



h 



Ox 1/T 



i\x-Xi\/u 
-i\x,\/u 



-i\x-xyu 



'^^ 1/T 

a) 



b) 



FIG. 24: Analytical structure of the integrand in Eq. (|E.2p 
in the complex plane of t[ and t\ . (a) The contour of integra- 
tion over t[ is closed upwards and runs around three branch 
cuts [Ji,2,3 are the corresponding terms in Eq. HE.4[) ]. (b) 
The contour of integration over ri for the first two terms in 
Eq. I|E.4|I is closed downwards and runs around two branch 
cuts. 



sinh(7rTz), J(a;, r, xi) is written as 



1/T 



cLti 



1/T 



si-°'= [i{T[ - ri)]si-"- [-i{T[ - Ti) 



where 



s{\xi\/u + iT'y)s{\xi\/u - ir{) 
s{\xi\/u + ZTi)s(|a;i|/u - in) 

s{\x — Xl\/u A- IT — ITl)'^ ""^^ 



Qe/2 



XI -1 



s{\x — Xi\/ U ^- IT — iT[) 

s{\x — Xl\/u — IT + ITl) 
s{\x — Xi\/u — IT + iT[) 



Qe/2 



X2-1 



(E.3) 



Xi = 

X2 = 
X 



s{xi/vp + iT[)s{xi/u + iTi) 

s{xi/u + iT[)s{xilvp + ITl) 
S[{X - Xl)/Vp + z(t - Tl)] 

s{{x — xi)/u + i(r — Tl)] 

s[{x — Xlj/u + i{T — t[)\ 
S[{X - Xl)/VF +i{T - t[)\ 



All functions B^^ in Eq. (jE.2p are understood as taken at 
zero disorder. The analytical structure of the integrand 
in the complex plane of ri and t[ is similar to that in the 
integrals analyzed in Sees. IVII ElVII PI Specifically, in 
the plane of t[ we have six vertical branch cuts starting at 
t[ = ±z|a;i|/M, t{ — zLi\x — xi\/u, and t{ = ti ±i/A, and 
going to Imr{ — ±zoo, Fig. [Mk . On top of that, there 
are singularities associated with the functions rj^ (j7.13p . 
Note that the strong singularity (a double pole if a — > 0) 
at t[ = Tl (if A is sent to infinity) that comes from the 

product of two free Green's functions g'f^ [ 0, ±(t{ — ti) ] 
is largely canceled by the vanishing of the rest of the in- 
tegrand at this same point. Evaluation of the integrals 
(|E.2[) at arbitrary a leads to rather cumbersome expres- 
sions which are not particularly interesting for our pur- 
poses. Let us extract the factor that renormalizes the 
impurity strength, 

X J{x,T,Xi) , 

and analyze the behavior of the factor J(a;,r, a;i) at 
small a. By introducing a short-hand notation s{z) = 



In Sec. IVII CI while calculating the disorder-induced 
correction to i?^^ to order 0(70;^), we neglected the dif- 
ference between vp and u in the interaction propagators. 
To compare with the result of Sec. IVII C) we neglect the 
difference between the two velocities here as well; in par- 
ticular, by putting xi = X2 = ^- We discuss the role of 
such factors in the calculation of the dephasing rate in 
Appendix |f1 

The integral over r{ in Eq. (|E.3[) can be represented as 
a sum of three terms, 

J = -^ I ' (Ji + J2 + J3) , (E.4) 
^ Jo 

each of which is the integral around a branch cut in the 
complex plane of t[, Fig. [Mb . Inspection of Eq. (|E.3p 
shows that J3, the term related to the branch cut start- 
ing at r{ = Tl -I- i/A, can be omitted for ae <C 1, since 
J3 ~ 0{al), which gives a subleading contribution to J 
[one power of ae comes from the difference of the inte- 
grand on two sides of the cut and two more come from 
the expansion of the expressions in the curly brackets in 
Eq. (jE.3p in r{ — ti]. The integrals Ji^2 are written for 
Qfe <C 1 as 
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Ji ~ a, 



nT 



dt[- 



,2(1 



X < 



/\x- 
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+ ir 
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\x—xi 



Qe/2 



In turn, each of the integrals of Ji and J2 over ri can be spht into a sum of two integrals along branch cuts in the 
complex plane of ri (Fig. [Mb), so that in total we have J ~ (Jia + Jn,) + (J2a + ^2b), where for ag ^ 1: 
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Jib — Q^e 



dt[ 



s ( i'l + 2 



ii + 2 



\x—xi I 



+ ir 



Uf 



\ ^ I'l^ I'F / \ J- Vf Vf / \ J- Vf Vf 



vf Vf 
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Vf 

Qe/2 
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J2a ~ a 



dti / dt[ -. 



irj 3 ^ 



ti + 



\xi\ . \x—xi I 



\x~xi\ _ \xi\ 
\ 1 Vf Vf 

^) 

./2 



e/2 



ti+2M),s(fJ,s(t;+2 



Ix— a;i I 



(E.7) 



r 



and J2b is obtained from Jia by interchanging a:i a: — 
xi. The terms Jiq and J2b on one hand and Ju and 
J2a on the other give essentially different contributions 
to J(a;,T, xi). The integrals Jia and J2&, as a function 
of the position of the impurity xi, are not exponentially 
small only in a close vicinity of width ~ u/T around 
xi = and xi = x, respectively. After averaging over xi, 
these "local" terms yield a contribution of order 0{al) 
which does not depend, in the leading approximation, on 
X and r. On the contrary, Jif, and J2a for small ae do not 



depend on 0:1 for xi between points xi = and 2:1 = x. 
The integration over xi for |a;|/?;F ^ ^/T yields then the 
leading contribution to Eq. (jE.ip . Only the latter terms 
are thus of importance to us. In the limit of small ae, 
their contribution to J is represented as 

2 



J{x,T,xi) ~aA — 



ttT 
2 



sinh 



dti 
"'0 



ttT Ui + t\ + -L-^ 

VF 



dt[ 

\x - Xi\ 
VF 



+ ir 
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+ (T^-r)}. (E.8) 

Analytically continuing J{x,t,xi) onto the imaginary 
axis of r, we observe that J{x,it,xi) is exponentially 
small for large positive values of (T/i'i;')(|a;i| + |a: — a;i|± 
vpt) and is a linear function of this sum when it is large 
and negative. Since we are interested in the limit of 
\x\/vf, \t\ ^ l/T", we can approximate J{x,it,Xi) as 

2 \ VF vp J 

X e - \x- xi\ +VFt) + {t^ -t) . (E.9) 

Integrating over xi we obtain 

/..,Z(.,..,.,) = -e-">""'-"M.,.<).(E.10) 

where b{x, it) is given by Eq. (|7.54p . We thus provided a 
direct verification of the RPA form of Eq. (|7.8[) by deriv- 
ing the first term in the expansion of exp(— _B_|__|_) in pow- 
ers of 7 without using the disorder-renormalized effective 
interaction (|7.20p . This procedure can be extended to 

higher powers of 7 and to the correlator 

An important point concerns the factors 
sinh^"°^^(7rTz), where z denotes the arguments of 
the functions s(z) in Eqs. (|E.6|) . (|E.7|) . These were 
omitted in the limit of small ae in Eq. ()E.8|) . However, 
even for ae <^ 1, this step cannot be justified for 
arbitrary x,xi, and r, e.g., if ae\x\T/vF ^ 1- What 
is important to us is that all these factors actually 
cancel each other after the analytical continuation onto 
the imaginary axis oi t it for large positive values 
of — |xi| — |a; — xi| ± ut. As Eq. (jE.9|) shows, only 
this range of t is relevant to J{x,it, xi), which is the 
quantity used in the calculation of the dephasing rate. 
The RPA for pairing of all four fields at points (0, 0) 
and (x, r), shown in Fig. I24[ namely of the combinations 
9+{x,t) - e^{x,T) and 6*+ (0,0) - 6'_(0,0) with each 
other, via the fermionic bubble with two backscatterings 
at X = xi can be reproduced in a similar way. 



APPENDIX F: CLASSICAL TRAJECTORY FOR 
THE COOPERON 

In this Appendix, we substantiate an important step 
in our derivation of the interaction-induced action Sc of 
the Cooperon within the functional-bosonization scheme. 
Specifically, we analyze the role of the factors '7±(a;,r), 
Eq. (|7.13|) . which contain both the Fermi velocity vp 
and the plasmon velocity u. We demonstrate that in 
the spinless case, considered in this paper, the velocity 
of electrons on the "Cooperon path" which contributes 
to ActwI is given by u and not hy vf- Put another way, 
the velocity of the interfering quasiclassical trajectories is 
renormalized by the interaction and coincides with that 
of the RPA propagators. As a result, one can omit the 



potentially dangerous factors r]±{x,T) in the calculation 
of the WL dephasing. 

The function C3 in Eq. (|7.6ip . which gives the three- 
impurity Cooperon loop, is explicitely written as 

C3(x,Tf - T,,Xi,X2,X3) = [-^— ) 

.l/T .l/T .l/T .l/T .l/T .l/T 

X dn dfi / dT2 / df2 / dra / df^ 
Jo Jo Jo Jo Jo Jo 

X <7i°^(x2 - Xi,T2 - Ti) e"'[«+(^2^^^)-''-(^^'^^)l 
X (X3 - X2,rs - T2) e^[«.(-3,r3)-.-(.3..3)] 
X - a;3,T-/ - T3) 

X gf{x2 - Xi,f2 - f^)e*[e+(:';2.r.)-6»_(:..,f2)] 

X gW(x3 - 2:2,7=3 - f2) e-['^+(--^3)-(^-(.3,^3)] 

X g^^\-X3,n-f3)) , (F.l) 

where we write the external times of the bubble as and 
Tf [in Eq. (|7.6ip . Tf — Ti = r] . To examine the role of the 
factors ?7±, it is convenient to introduce the function 

Jo Jo 

by integrating first over the external times. 

The averaging (. . .) over the fluctuating fields 0±{x, t) 
yields the factor 

exp[-2(M//-M/b)] , 

defined in Eqs. ([7\37l) - l[7?39ll . in the integrand of 
Eq. (jF.ip . In the absence of interaction {Mff = Mff, = 
0), the integrals over the internal times are determined 
by the poles of the free Green's functions [specified by 
Eq. (IF.4|) below with u — > vf]- For given coordinates xi, 
X2, X3, and X, we refer to this set of times as the "classical 
trajectory" of the Cooperon with velocity vp- 

Each interaction-induced term M{N,N') in 
Eqs. (|7.37p . (|7.38p contains the combination + 
and hence depends on both u and vp. If one neglects 
disorder in the RPA propagators ("clean RPA"), the 
factor exp[~M{N, N')] is represented as 

^-M(N,N') _ ^-a^C(xN~x^i .TN~Tf^,) 

X C{xn - Xn',Tn ~ Tn') , (F.2) 

where Cix,T) — r]+{x,T)ri-{x,T), ae = 1 — K, 
and the functions C{x,t) and 'ri±[x,T) are given by 
Eqs. (|7.12p . (|7.13p . The dependence 0x1 vp comes from 
the factors rj±. The "dirty RPA" does not affect the 
latter. We thus see that each pair of impurities N and 
N' in Eq. (jF.ip generates two poles associated with the 
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functions r]±: 



(?) 

sinh [ t:T{x/vf + it) ] sinh [ ■kT{x /vp — it) 
sinh ['kT{x/u + ir) ] sinh ['kT{x/u — ir) ] 



The singularities in Eq. (jF.ip coming from the functions 
C,{x^t) are analogous to those in the bare Green's func- 
tions g''^{x,T). However, they can appear either as "m- 
poles" at IT — ±x/u or as "tiF-poles" at it = ±x/vfj de- 
pending on the sign of the factor (—1)^+^ or (—1)^+^ 
in Eqs. (ESZD, (17351). 

We can readily integrate over and Tf the free Green's 
functions that connect to the external vertices. Since 
there are no interaction-induced factors at the external 
vertices, the integration yields again the free Green's 
functions. The function C3 then becomes 



V 2 / 



l/T 

(0)/ 
(0)/ 

gL\xi 



l/T 

dfi 



l/T 
dT2 



l/T 

df2 



Jo Jo 

X3,Tl ~f3)C{xi 



(0)/ 

(0)/ 
9-- [X3 



C\X2 



Xl,T2 

X2,T3 

X3,Tl 

Xl,f2 

X2,f3 

X3,Tl 
X2,Ti 
Xl,T2 



Tl) C{X2 
T2) C{X3 

T3) C{xi 

n)C{x2 

T2) C{X3 



10 
X3,Ti 

Xl,T2 

X2,T3 

X3,Tl 

Xi,f2 



l/T pl/T 

dT3 / dfs 
Jo 



r3) 
ri) 
r2) 
r3) 
fi) 
r2) 



X2,T3 

T3)C^{X1 -X3,fi -fs) 
T2) C^{X3 - X2,T3 - f2) 
n) C"^(a;2 - X3,T2 - fa) 
C(0,Ti~fi)C(0,T2-f2)C(0,r3-f3) 

2:2, Xs, Tl, T2, Ta, fi, f2, fa) 
W{xi,Ti,X3,f3,nn) 

iy(a;-xa,T3,x-a;i,fi,-f]„) , (F.3) 



where C(0, r) = (vp/u)'^, the function R combines all the 
factors coming from the first factor in Eq. (jF.2[) . and 



W{xi,Ti,X3,T3,nn) 

= e 
+ e 



\flnXi 



'''-^'[Q{nn)Q(xi)-e{-nn)Q{-xi)] 
^""^MQiaoeia^s) - e(~r!„)0(-2:3)] 



In the single-particle Green's function, the Fermi 
velocity vp drops out completely in the product 
(7±^(a;,r)77±(x,T). One can see that each free Green's 



function in Eq. (jF.3[) is multiplied by a corresponding 
factor ri±, so that all the tij^-poles associated with the 
free Green's functions in Eq. (jF.3|) are replaced by u- 
poles. One might think that vp disappears completely 



from the problem, and hence the integration over the in- 
ternal times in Eq. (|F.3p straightforwardly yields a new 
classical trajectory defined by 

iri = —xi/u , 

iT2 = {x2 - 2xi)/u , 

ITS = (2x2 -2xi - X3)/u , 

if 3 = -X3/U , 

if2 = {X2 - 2X3) /u , 

ifi = {2x2 - 2x3 ~ xi) / u . (F.4) 

However, the impurity vertices are all interconnected by 
the interaction propagators. Among the latter, there 
are propagators that connect vertices which are not 
connected by the single-particle Green's functions [e.g., 
points (xi,Ti) and (a;2,f2)]. As a result, the w^-poles in 
the corresponding functions ri± remain uncompensated. 

In the noninteracting case, the sequence of points in 
real space 



1 



1 



1 



(F.5) 



forms a closed loop with segments given by the Green's 
functions g^^*^ . Now that we have the interaction-induced 
functions ?7± that connect every pair of impurities, other 
trajectories can contribute to C3, whose segments are 
given either by the Green's function or by the function 
•q±. In view of this, the notion of a single classical tra- 
jectory contributing to Acr„i requires a more accurate 
justification. As we show below, the interaction itself 
"chooses" a unique classical trajectory which is defined 
in Eq. (|F.4p . i.e., the one that is determined by a com- 
bination of the Green's functions with the renormalized 
velocity vp — > u. The contributions of all other trajec- 
tories are exponentially suppressed. Note that within 
the fermionic path-integral approach used in Sec. lVIi the 
"correct" velocity is chosen automatically on the saddle- 
point trajectory of the total action. 

When analyzing the role of the factors •q±, we can, 
in the first approximation, neglect other interaction- 
induced factors in Eq. (jF.ip by setting the exponent 
tte = in Eq. ([Q]) . i.e., by putting i? = 1 in Eq. ([Q]) 
while keeping the difference between vp and u in the 
functions C- This simplifies the analytical structure of 
the integrand in Eq. (|F.3p by reducing the integration 
along branch cuts to the integration around simple poles 
only. The integration around all the poles then yields C3 
represented as a sum over a whole set of trajectories, not 
only the one given by Eq. (jF.Sp . 

The integrand of Eq. (jF.3p contains products of the 

type g'£\x,it)C,{x,it) which, as a function of t, have 



poles on the classical trajectory (jF.4p at t 
a residue equal to 

(0). 



±x/u with 



[ut ^ {x + u sgnz/A) ] 9± {x, it) ({x, it)\t-t 



^x/u 



1 Vp 

— — exp 

Zm u 



Vp 



UVp 



ttT\x\ 



(F.6) 
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for \x\/u ^ 1/T. Writing the remaining factors C which 
are not singular at this pole as 

C(x,^^) ^ (!^)'exp[7rT(|t;| + |r|-|t!t|-|i^|)] 
X sgnt'; sgnT sgni!^ sgnT (F.7) 

for |t!^'"| > 1/T, where t"^^t± x/v and f^=t± x/u, 
we observe that all the exponential factors coming from 
Eqs. (|F.6p . (IF.7P cancel out on the classical trajectory 

dEl- 

On the contrary, any other sequence of poles has at 
least some of the segments of the trajectory passed with 
velocity vp and no such compensation occurs, which 
leads to an exponential decay of the contribution to 
C3. The exponent depends on a particular sequence of 
poles but for any sequence is at least of the order of 



{\u — vp\/vp)Ttc. One sees that the contribution of any 
of the "nonclassical" trajectories is exponentially sup- 
pressed either by the "Golden-Rule dephasing" or sim- 
ply by the thermal smearing. The role of the factors ri± 
is thus two-fold: they renormalize the velocity on the 
"classical" Cooperon trajectory {vp — > u) and suppress 
the contribution of all nonclassical trajectories. Taking 
into account the factor R in Eq. (jP.Sp within the "clean 
RPA" does not give any decay on this classical trajec- 
tory either. Including disorder within the "dirty RPA" 
(Sec. IVII CI and Appendix [E} in the factor R yields an 
exponential decay on the classical trajectory as well, as 
discussed in Sec. IVIIFl but this decay is much weaker 
than on the nonclassical trajectories, so that the latter 
can indeed be neglected in the calculation of ActwI. 
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